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MAGNETO-OPTIC  MATERIALS  FOR  BIASING 
RING  LASER  GYROS  -  REPORT  NO.  3 
(COMPUTER  MODEL  FOR  EVALUATING  SCATTERING 
FROM  MULTI-LAYER  DIELECTRIC  THIN  FILM  STRUCTURES 
CONTAINING  A  MAGNETIC  LAYER) 

I.  INTRODUCTION  AND  OVERVIEW 

The  counterpropagating  beams  in  a  ring  laser  can  be  employed  to  sense  physical  rotation  and 
hence  form  the  basis  for  a  ring  laser  gyroscope  (RLG)  [1].  Via  the  Sagnac  effect,  the  two  similar  coun¬ 
terpropagating  beams  oscillate  at  slightly  different  frequencies  with  the  difference  frequency  given 
by 

Ai»-<Mn/LA  (1) 

where  A  is  the  enclosed  area  and  L  the  optical  length  of  the  RLG,  O  is  the  rotation  rate  and  A,  the 
laser  optical  wavelength.  When  0  becomes  too  small  the  two  beams  can  lock  at  a  common  frequency 
with  consequent  loss  of  rate  information.  This  lock-in  problem  arises  from  optical  backscatter  which 
couples  the  beams.  As  a  result,  it  is  important  to  keep  the  optical  path  as  free  of  additional  elements  as 
possible. 

Because  of  their  non-reciprocal  propagation  properties,  magneto-optical  ( M-0 )  elements  can 
introduce  a  controllable  frequency  difference  or  magnetic  bias  and  thus  avoid  the  lock-in  problem 
in  an  elegant  manner.  The  magneto-optical  effects  which  are  linear  in  the  magnetization  of  the  material 
and  hence  are  useful  for  this  application  are:  (1)  the  Faraday  effect  in  transmission  and  (2)  the  mag¬ 
netic  Kerr  effect  in  reflection.  The  Kerr  effect  has  three  subcases:  transverse,  polar  and  longitudinal, 
depending  on  how  the  magnetization  is  aligned  with  respect  to  the  plane  of  the  sample  and  the  plane  of 
incidence  of  the  light  beams.  These  M—  0  effects  are  discussed  in  Report  No.  1  [2]  of  this  series  where 
it  is  shown  that  they  each  introduce  a  differential  phase  shift  A4>  between  properly  chosen,  similar 

counterpropagating  beams.  The  magnetic  bias  induced  is  given  by  [1] 
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Av4  »  (A4>/2 ir)c/L  (2) 

where  A4>  is  in  radians  and  c  is  the  velocity  of  light.  A  review  article  on  the  applications  of  magneto¬ 
optics  in  ring  laser  gyroscopes  has  recently  been  prepared  [3]  and  provides  a  concise  unified  treatment 
of  the  subject. 

It  is  advantageous  to  use  the  Kerr  M-0  effects  for  achieving  magnetic  bias  because  then  the 
M—O  element  can  be  used  as  one  of  the  mirrors  which  define  the  RLG  cavity.  This  means  that  it  is 
not  necessary  to  place  any  additional  elements  in  the  beam  as  would  be  the  case  if  one  used  the  Faraday 
effect.  As  a  result  the  problems  related  with  scattering  from  extra  surfaces  are  eliminated. 

The  reflectivity  of  typical  magnetic  metal  films  is  only  in  the  range  0.4-0.7  for  the  He-Ne 
wavelengths  0.63  ^m  and  1.15  /xm  used  in  RLGs.  These  low  reflectivities  are  inconsistent  with  the 
optical  loss  requirements  of  RLGs  (<5%  at  1.15  pm  and  <  1-2%  at  0.63  /j.m)  and  with  the  capabilities 
of  multilayer  dielectric  (MLD)  mirrors  which  have  total  optical  losses  of  only  0.02%  for  laser  quality 
mirrors.  In  order  to  reduce  the  loss  associated  with  the  M-0  mirror  one  can  overcoat  the  opaque 
metallic  magnetic  layer  with  an  MLD  stack  which  is  tuned  to  be  highly  reflecting  at  the  wavelength  and 
angle  of  incidence  used.  Similarly,  if  a  transparent  magnetic  material  is  used,  it  can  be  incorporated 
into  an  MLD  stack,  again  to  achieve  adequate  reflectivity.  RLG  mirrors  based  on  both  opaque  metallic 
and  transparent  garnet  magnetic  films  have  been  discussed  by  workers  at  Sperry  [4.5]. 

The  above  considerations  indicate  that  it  is  highly  desirable  to  have  a  unified  method  of  evaluating 
the  important  optical  and  magneto-optical  parameters  which  characterize  a  given  multilayer  magnetic- 
dielectric  thin  film  composite.  A  computer  model  which  permits  various  RLG  M-0  multilayer  mirrors 
to  be  evaluated  and  compared  quantitatively  would  be  a  very  useful  design  tool.  With  this  goal  in  mind 
we  have  produced  such  a  computer  model  whose  characteristics,  structure,  and  details  are  set  forth  in 
the  bulk  of  this  report. 
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In  Chapter  II  the  mathematical  basis  of  the  method,  which  is  the  scattering  matrix  approach 
developed  by  Hunt  [6]  for  multilayer  optical  structures  containing  a  single  magnetic  layer,  is  described. 
Chapter  III  shows  how  these  ideas  have  been  incorporated  into  a  computer  model  (program)  which  can 
calculate  the  differential  phase  shift,  the  total  reflectivity,  the  differential  reflectivity  (for  similar  coun- 
terpropagating  beams),  and  the  mode  mixing,  if  any,  for  a  general  multilayer  structure.  All  three  Kerr 
effects  (transverse,  polar  and  longitudinal)  can  be  handled  easily  by  the  model  which  also  takes  into 
account  Faraday  rotation  in  the  magnetic  layer.  In  Chapter  IV,  results  from  calculations  on  several  typ¬ 
ical  structures  which  might  be  used  as  magnetic  bias  mirrors  in  RLG  applications  are  illustrated. 
Finally,  the  Fortran  listings  of  the  main  program  and  its  subroutines  are  given  in  the  Appendix. 

The  calculations  given  in  Chapter  IV  indicate  the  superiority  of  mirrors  based  on  magnetic  garnets 
to  those  based  on  transition  metals  such  as  iron.  They  also  show,  however,  that  MnBi  polar  Kerr  mir¬ 
rors  are  a  viable  alternative  to  the  garnet  mirrors.  It  appears  likely  that  the  choice  will  depend  on  the 
cost  and  reproducibility  factors  in  the  fabrication  of  garnet  or  MnBi  mirrors. 

Readers  may  obtain  Fortran  versions  of  the  computer  program  on  their  flexible  disc  media  by  con¬ 
tacting  one  of  the  authors  at  the  Naval  Research  Laboratory,  phone  number  (202)-767-3603. 
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II.  MATHEMATICAL  BASIS  FOR  THE  MODEL 

A  general  method  of  calculation  is  required  if  one  is  to  handle  a  wide  variety  of  multilayer  dielec¬ 
tric  structures  which  contain  a  layer  having  magneto-optical  properties.  Fortunately,  general  scattering 
matrix  methods  for  evaluating  the  magneto-optic  scattering  from  such  structures  were  developed  in  the 
late  1960’s  by  Smith  [7]  and  Hunt.  [6]  Our  embodiment  of  the  scattering  matrix  techniques  into  a  com¬ 
puter  program  possessing  a  high  degree  of  generality  is  based  on  Hunt’s  methods.  In  this  chapter  we 
first  outline  the  method  in  enough  detail  to  make  our  computer  program  understandable  and  note  an 
error  in  Hunt’s  paper  which  is  corrected  in  the  program  used  here.  That  error,  in  the  transverse  Ken 
Fresnel  coefficients  and  t„,  was  pointed  out  to  us  by  R.E.  McClure  and  the  correct  derivation  can  be 
found  in  his  recent  paper  [8].  Finally,  we  write  out  explicit  expressions  for  the  differential  phase  shift, 
reflectivity,  differential  reflectivity,  and  mode  mixing  for  ring  laser  gyro  M—O  minors  using  each  of  the 
three  Ken  effects. 

In  Fig.  1  a  general  thin  film  structure  is  shown.  It  consists  of  L  layers  each  having  a  specified 
thickness  dk  and  index  of  refraction  nk.  There  is  a  single  magnetic  layer  characterized  by  du,  nM  and 
the  magneto-optical  coefficient  Q.  Note  that  and  Q  are  complex  quantities  while  the  other  parame¬ 
ters  are  real. 

Let  Ef  •*  be  the  total  electric  vector  of  a  right  propagating  p- polarized  beam  at  the  k  interface,  in 
the  k  -  1  medium.  Similarly,  £,“  *  is  the  total  electric  vector  of  a  left  propagating  s-polarized  beam  at 
the  k  interface  in  the  k  -  1  medium.  Note  that  the  k  interface  separates  medium  k  -  1  from  medium 
k.  Then  the  4  x  4  scattering  matrix  (Sk)  which  describes  the  effect  of  interface  k  and  propagation  in 
the  k  medium  obeys  the  following  equation  (6], 


£-•* 

£„-*+l 

E,-k 

£f*+l 

V’* 

-  (S*) 

£„+*+1 

£,*•* 

£  +  .*  +  1 
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Hunt  (6]  shows  if  the  k  medium  has  magneto-optical  properties  then  the  front-surface  scattering 
matrix  {S’*)  is  given  by 


(S*0  - 


[?  -  /•*(/*)- *7*]  Tk  r*(r*)~*  Tkl 

—uk)~lTkTk  (r*)~l  rr1 


(4) 


where  r\  7*.  r  and  r  are  2  x  2  matrices  of  Fresnel  coefficients.  The  Fresnel  matrices  are  defined  by 


r* 

rPP 

rk 

’pi 

Tk 

rpp 

7k 

rpi 

rk 

r*P 

Tk 

'is 

.  r- 

7k 

rsp 

7k 

'ss 

(5) 


where  the  Fresnel  reflection  coefficients  at  the  k  interface  are  rjj  for  an  incident  beam  propagating  in  the 
dielectric  toward  the  magnetic  medium  and  are  Tj*  for  an  incident  beam  propagating  in  the  magnetic 
medium  toward  the  dielectric.  Note  that  p  and  s  refer  to  the  plane  of  polarization  [see  Report  No.  1, 
Fig.  2  and  Eq.  (11-6)1  [2].  The  tk  and  7*  matrices  contain  the  corresponding  transmission  Fresnel 
coefficients. 


In  Eq.  (2)  Tk  is  a  propagation  matrix  which  describes  the  normal  phase  retardation  and  the  Fara¬ 
day  effect  for  a  wave  (£„*,  E*)  propagating  through  medium  k  of  thickness  dk.  Hunt,  based  on  the 
work  of  Robinson,  [9]  shows  that  Tk  has  the  form  [6] 


where 


Tk 


cos</>* 

-  sin  t l>k 


sin  t/>* 
cos  l llk 


(6) 


<6*  —  2irnkdk  cos  0*/k  (7) 

and  t/i*  depends  on  the  magneto-optical  effect  (transverse,  longitudinal  or  polar)  with  which  we  are 
dealing  and  can  be  found  in  Table  II.  The  angle  between  the  beam  direction  in  medium  k  and  the  frlh 
interface  normal  is  0*.  Note  that  Tk  is  the  same  for  waves  propagating  in  the  right  and  left  directions. 
When  the  last  layer  L  is  opaque,  or  is  of  indefinite  (but  large)  thickness,  or  if  propagation  effect  for  this 
layer  are  otherwise  not  important,  one  can  set  dL  *■  0.  This  converts  TL  into  a  unit  matrix  and  elim¬ 
inates  the  propagation  effects  while  keeping  the  interface  effects.  Since  it*  —  nm  is  complex,  so  is  <6* 
(plus  0*  and  t //*).  This  takes  into  account  optical  absorption  in  the  magnetic  medium. 
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The  Fresnel  coefficients  for  a  M—O  material  in  each  of  the  M—0  modes  are  derived  by  Hunt  and 
we  give  these  expressions  for  rfj,  rjj,  tfi  and  and  <fik  in  Tables  I  and  II.  As  can  be  seen,  the  Fresnel 
coefficients  for  the  k  interface  involve  nk-u  nk,  0*-i<  0k  and  Q.  We  point  out  that  all  off-diagonal  com¬ 
ponents  of  r*.  7*,  r*  and  7*  are  proportional  to  Q  (or  all  vanish  for  the  transverse  mode).  Since  |(?| 
<  <  1,  one  can  neglect  second  order  terms  in  Q.  This  means,  for  example,  that 


(/)-' 


tpp  tpt  1/ tpp  tpj  tppta 

t sp  t$s  “(jp/  tpp^u  1/ 4j 


where  we  have  suppressed  the  superscript  k. 


(8) 


Table  I.  Single-Interface  Diagonal  Fresnel  Coefficients 


Element 

Longitudinal 

and  Polar 

Transverse 

n2Ci  -  rt\c2 

n2cx(l  +  iQs-t/cJ  -  n\C2 

rpp 

n2ci  +  nxc2 

n2c  i(l  +  iQs-JcJ  +  nxc2 

ntCi  -  n2c2 

»ict  -  n2c2 

rss 

n,ci  +  n2c2 

«ic,  +  rt2c2 

tpp 

W«2>  (1  +  rw) 

(n\/n2)(l  +  r„) 

ta 

1  +  ra 

1  +  rra 

n2C\  -  nic2(l  +  iQsi/cJ 

fpp 

~rpp 

n2c\  +  rtic2(l  -  /Qs2/c2) 

~ra 

tpp 

(«2/«t)(l  -  r„) 

(n^nx)  (1  +  ?„) 

1  -  ra 

1  -  r*s 

<t> 

ford^c-JX 

ft),  n 2  —  refractive  indices  of  the  dielectric  and  magnetic  media,  respec¬ 
tively. 

C|,  c j  —  direction  cosines  for  the  angle  of  refraction  in  the  dielectric  and 
magnetic  media,  respectively. 

s 2  —  sine  of  the  refraction  angle  in  the  magnetic  medium. 
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TABLE  II.  Single-Interface  Off-Diagonal  Fresnel  Coefficients 


Element 

Polar 

Longitudinal 

rv 

j  Qtppt^sJUxS-i) 

J  QtpptxSx/iCxCi) 

rps 

~rsp 

rv 

(»• 

{njn2)rps(  1  -  2clc2/tpp) 

(ni/n2)rps(l  +  2c\C2/tpp) 

7y> 

~rpprsplfpp 

~rppripl  (pp 

tpsl  *33 

rsi  tps/  ta 

~rpprsp/tpp 

~rpprJ  tpp 

<p. 

~rssrpJ  *a 

~rurps/  la 

* 

ird2n2Q/\ 

irdj^Qisi/  c2)/k 

«i,  —  refractive  indices  of  the  dielectric  and  magnetic  media  respec¬ 

tively. 

C|,  c2(j,.  Jj)  —  cosine  (sine)  of  the  refraction  angle  in  the  dielectric  and 
magnetic  media,  respectively. 

One  can  now  write  out  explicit  expressions  for  (S^poian  (•S‘M)|0ng  and  (S*0trin ,  and  these  are  given 
in  Tables  III  and  IV.  Note  that  and  (Sw),ont  have  the  same  form  while  (5*0trins  is  simpler 

because  of  the  lack  of  mode  mixing  in  the  transverse  mode.  In  Table  III,  terms  of  the  form  Q  sin  4>k 
have  been  neglected  as  being  of  order  Q1.  Although  this  approximation  could  be  invalid  for  a  thick 
transparent  M-0  slab,  we  anticipate  that  *k  <  <  1  for  practical  ring  laser  gyro  applications  and  hence 
the  approximation  is  valid. 


KREBS  AND  MAISCH 


TABLE  III.  Front  Surface  Magneto-Optical  Scattering  Matrix  —  Polar  and  Longitudinal  Kerr  Effects 


*dt„ 

*s/t„ 

[(fpp^pj  ips  f pp ) C j  { tpp  )  FppSl 

Qc/ta 

(rvc  +  rat„s)IUppta<^) 

rac/  (fu4>) 

<Pr„c/t„ 

QTppSjtpp 

cKw D) 

-(/pjC  +  tas)/ (tpptxQ) 

-Qr^s/u j 

♦  r^c/fts 

(fpp-T  "  tspC}/ 

c/  (/»<!>) 

—  e'*.  c  —  cos  <b.  s  —  sin  < b 

The  superscript  k(—M)  has  been  dropped  on  the  Fresnel  coefficients  as  well  as  on  <6  and  ib. 


TABLE  IV.  Front  Surface  Magneto-Optical  Scattering  Matrix  —  Transverse  Kerr  Effect 


(fppfpp  ~  rwrpp)tl >/fw 

0 

r pp/  ( fpp*P  ) 

0 

0 

0 

rj  (f„<P) 

(5")- 

—7w<t>/fpp 

0 

I/Opp*) 

0 

0 

0 

<t>  -  e‘ * 

The  superscript  k(—M)  has  been  dropped  on  the  Fresnel  coefficient  as  well  as  on  <b. 


In  order  to  obtain  the  back-surface  scattering  matrix  (Sv+1),  one  must  reverse  the  roles  of  rk  and 
r*  with  respect  to  7*  and  7*.  Furthermore,  because  of  the  orientation  of  the  magnetization  vector  with 
respect  to  the  front  and  back  surface  normals,  we  must  make  the  substitutions 

/■((?) -7(-Q),  r(Q)-*r(-Q) 

r((2)-7(-0),  l(Q)-ti-Q) 

for  the  polar  and  transverse  cases.  For  the  longitudinal  case  this  sign  reversal  does  not  occur  and  we 


r(Q)  —  7(Q),  etc. 

Upon  making  these  substitutions  in  Tables  III  and  IV  we  get  the  back-surface  results  shown  in  Tables  V 
and  VI. 
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TABLE  V.  Back  Surface  Magneto-Optical  Scattering  Matrix  —  Polar  and  Longitudinal  Kerr  Effects 


(Sw+1) 


*pp 

-<J*r*/  (App/^) 

* rpp/lpp 

<7*^/ (A,/*) 


G<t>rpp/  (A  tpp) 

*/7„ 

G<t>rps/  (A„7W) 
*  r */  r* 


fpp!  ( ipp* ) 

G  fppT  */  (A pp  /** ) 

l/(7wd>) 

G  Ppp  P^p/  (A pp  r** ) 


G  rssr  pp!  (A*/^*) 

G  F’ssl’ps/  (A*rpp*) 

I/O**) 


4>  - 

The  superscript  Id—M+l)  has  been  dropped  on  the  Fresnel  coefficients  and  <t>. 
Aa  -  (7*)_2.  1-(7„)J  _ 

('pptps- W,,) 

G  —  +  1  (Polar),  G  —  —  1  (Longitudinal) 


TABLE  VI.  Back  Surface  Magneto-Optical  Scattering  Matrix  —  Transverse  Kerr  Effect 


(S‘M+1) 


[U-0)/„p(-O>  -  r.A-Q^A-QWJj-Q)  0  7pA-Q)/[tpp{-Q)^}  0 


-rpp(-Q)<t>/7pp(-Q) 

0 


*/f* 

0 


i/[fpp(-0)*l 

o 


^/Oss*) 

0 

1/(7**) 


<t>  -  e'<* 

The  superscript  k  (—  Af+1)  has  been  dropped  on  the  Fresnel  coefficients  and  <t>. 

Finally,  we  note  that  if  we  are  dealing  with  an  interface  at  which  both  media  are  nonmagnetic 
dielectrics,  then  Eq.  (S)  reduces  to  the  simple  form 


($*)*,.  - 


**/£ 

0 

r*,/(r***) 

0 

0 

**/4 

0 

4/(4**) 

**r* /tk 

^  rppi  ‘pp 

0 

1/(4**) 

0 

0 

**4/  4 

0 

1/(4**) 

(9) 


where  **  —  e'*k  and  the  Fresnel  coefficients  can  be  obtained  from  Table  I  with  Q  —  0. 
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The  overall  scattering  matrix  for  the  general  multilayer  magnetic-dielectric  thin-film  structure 
shown  in  Fig.  1  is  then 

(s)  -  (s‘)(s2)  -  -  (SL~l)(SL)  (10) 

in  which  all  but  (5**)  and  (SM+{)  are  dielectric-dielectric  matrices  of  the  form  Eq.  (7).  Thus  we  have 


(11) 


The  scattering  matrix  (5)  contains  all  the  desired  information  about  the  optical  and  magneto-optical 
properties  of  the  multilayer  structure.  One  can  obtain  the  Fresnel  reflection  and  transmission 
coefficients  for  the  entire  stack  as  well  as  the  Kerr  rotations  K„  and  Ks  and  Faraday  rotations  Fp  and  Fs. 
(All  of  these  quantities,  like  the  components  of  (5),  are  in  general  complex  numbers.)  The  equations 
which  relate  these  stack  Fresnel  and  rotation  parameters  to  the  components  of  (S)  are  derived  by  Hunt 
[61  and  are  summarized  in  Table  VII. 


Ep-° 

E,'L 

E,-° 

e:l 

Epa 

-  (5) 

e;l 

£,+,° 

E?L 

TABLE  VII.  Fresnel  Coefficients  for  Entire  Stack 


'pp 


”  (S«S  13  —  S^S^)/  (S32S44  —  S34S43) 
“  (S23S34  —  524S33)/  (534S43  —  S33S44) 
—  (5345|3  -  S14S33)/  (S34S43  —  S33544) 
“  (S23S44  —  S24543)/  (S33S44  —  S34S43) 


tan  K„  -  Rsp/Rpp  tan  K,  -  R^/R^ 


1pp 


‘St 

tan  F„ 


■S44/  (‘S33S44  “  S34S43) 

S33/  (533544  —  S34S43)  ■ 
—S3J  (S33S44  —  S34S43 ) 
—S43/  ( S33S44  —  534543) 

p!  Fpp  tan  Fs  ™  T^j  Ta 
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The  ring  laser  gyro  supports  counterpropagating  models  of  the  same  polarization (s).  The  desired 
optical  and  magneto-optical  characteristics  of  a  magnetic-dielectric  multilayer  structure  are:  1)  the 
differential  phase  shift  A4>  between  similar  counterpropagating  beams,  2)  the  reflectivity  R  of  the  stack 
for  these  beams,  3)  the  differential  reflectivity  AR,  4)  the  figure  of  merit  FM  =  A«fr/(1  -  R),  and  S) 
the  mode  mixing  e  =  ratio  of  the  improper  to  the  desired  polarization  amplitudes.  Since  the  p  and  s 
polarizations  are  eigenmodes  in  the  transverse  case,  «  =  0  for  transverse  M—O  Kerr  mirrors.  We  now 
consider  each  of  the  A/—  O  cases  in  turn. 

A.  Transverse  Kerr  Effect 

1)  Differential  Phase  shift 

A <t>T  SArg[E~(Q))  -  Arg[E~(-Q )) 

-  Arg{R„,(Q))  -  Arg[Rpp(-Q))  (12) 

where  Arg[X)  means  the  argument  of  the  complex  quantity  X  and  the  R,j  are  the  reflection  Fresnel 

coefficients  for  the  entire  stack  (obtained  from  the  Su  via  Table  VII).  We  have  used  the  fact  that  the 

magnetization  (and  Q)  looks  reversed  to  the  two  counterpropagating  ^-polarized  beams. 

2)  Reflectivity 

lff„(0)!2  -  y{ltf„(G)l2+  l*w(-0)l2}  (13) 

3)  Differential  Reflectivity 

A Rr  =  (|R„«?)IJ  -  |/?,,(-Q)l2/l/?,,(0)|J.  (14) 

The  calculation  of  FMT  is  obvious  and  c  T  =  0. 

B.  Polar  Kerr  Effect 

It  is  more  natural  to  deal  with  circularly  polarized  beams  in  this  case  so  we  use  the  circular 
Fresnel  coefficients  for  the  suck,  viz.  Rrk,  Rll,  Rlk  and  Rrl.  It  is  easy  to  show  that  these  can  be 
expressed  in  terms  of  the  R etc.  as  follows. 
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Rrr  Rrl 

_  1 

Rlr  Rll 

“  2 

(15) 


(R„  +  Ra  -  iRf,  +  iRv)  (R„  -  Ru  +  &„+  iR„) 

(R„  -  Ra-  iR,,-  iR„ )  (R„  +  R„  +  iR„  -  iRv) 

Since  Ra  =*  - R „  and  R„  =*  R„,  we  see  that  \Rrr\,  \RU\  <<  lR*zl,  iRtal,  ie.,  for  the  most  part 

R  —  L  and  L  — ►  R  upon  reflection.  However,  since  Rrr  and  RLL  *  0,  there  is  mode  mixing  and 
induced  ellipticity. 


1)  Differential  Phase  shift 

A<t>,  5  ArglER)L.m  -  ArglEC)^ 

-  Arg[RRL)  -  Arg[RLR)  (16) 

2)  Reflectivity 

\Rrl(0)\2-  |Rl*(0)|2 

-  |  !*„«»  -  R„(0)|2  (17) 

3)  Differential  Reflectivity 

AR,  -  URi*]2  -  !R*J2)/|R*£(0)|2  (18) 

4)  Mode  Mixing 

*PR  “  I  Rrr/ &LR  I 
*«.  “  \Rll/Rrl\. 

A  ring  laser  gyro  using  a  polar  Kerr  effect  M-  O  mirror  must  have  an  even  number  of  mirrors  in 
the  optical  path.  It  then  will  support  two  pairs  of  counterpropagating  beams,  one  pair  right  and  one  pair 

left  circularly  polarized  (as  referenced  to  a  fixed  point  in  one  arm  of  the  laser).  The  differential  phase 
shifts  for  the  RCP  and  LCP  pairs  are  equal  and  opposite. 

C.  Longitudinal  Kerr  Effect 

We  have  carried  out  a  Jones  matrix  analysis  of  longitudinal  M—O  Kerr  effect  mirrors  in  a  ring 
laser  gyro  system.  This  analysis  indicates  that  for  a  three  (or  other  odd  number)  mirror  laser  the 


(19a) 

(19b) 
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differential  phase  shift  would  depend  on  Q 2  rather  than  Q  so  that  such  a  case  would  not  be  useful. 
However,  in  a  four  (or  other  even)  mirror  system  there  are  differential  phase  shifts  linearly  dependent 
on  Q.  The  approximate  normal  modes  of  such  a  system  are  orthogonal  linearly  polarized  beams  of  the 
form 


Ea  "  Ep  +  Ef  ED- Ep-  E, 

where  the  E,  is  the  electric  vector.  The  longitudinal  analysis  below  is  therefore  cast  in  terms  of  EA  and 
Ed.  The  Fresnel  coefficients  are  then 


Raa  Rad 

1 

(Rpp  +  Ru  +  Rps  +  Rsp)  ( Rpp  -  Rss~  Rps  +  *„> 

Rda  Rdd 

"  2 

(Rpp  R* s  Rps  *”  Rsp)  (Rpp  Rss  Rps  Rsp) 

Again,  since  /?„=*-  Rpp  and  Rps  =  -  Rv,  we  get  |J?4j>|,  \Rda  I  >>  l^ool  30(1  hence  the 

polarization  essentially  switches  A  — *  D  and  D  -*  A  at  each  mirror. 


1)  Differential  Phase  Shift 

~  Arg[RD4 )  (21) 

2)  Reflectivity 

|/t.<D(0)|2-  |JtM(0)i2 

-  i|Rw(0)  - /?o(0)|J  (22) 

3)  Differential  Reflectivity 

LRl  “  (l*4ol2  -  l/?D4l2)/l^p(0)|2  (23) 

4)  Mode  Mixing 

eLA  *  \Raa/Rda\  (24a) 

*  ld  m  \Rdd/Rad\ •  (24b) 

The  computer  program  is  set  up  to  actually  carry  out  the  calculations  sketched  above,  keeping 
track  of  which  mode  (transverse,  polar  or  longitudinal)  is  being  used  and  hence  selecting  the  correct 
equations  to  evaluate  the  desired  quantities. 
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III.  THE  COMPUTER  MODEL 

The  basic  calculation  which  the  computer  model  must  make  is  to  determine  the  4  x  4  scattering 
matrix  ( S )  for  a  multilayer  magnetic-dielectric  structure  of  the  type  shown  in  Fig.  1.  As  shown  by  Eq. 
(10)  in  the  last  chapter,  this  means  one  must  evaluate  the  scattering  matrices  (5*)  for  each  interface- 
layer  and  perform  the  indicated  matrix  mutliplication.  Based  on  (5)  one  can  calculate  Rpp,  R„,  Rv  and 
Ru  for  the  entire  M-0  structure.  These  are  sufficient  to  evaluate  the  differential  phase  shift, 
reflectivity,  differential  reflectivity  and  mode  mixing,  as  well  as  the  figure  of  merit,  for  any  of  the  three 
M—O  modes  which  might  be  used  in  a  ring  laser  gyro  magneto-optical  mirror. 

The  simplest  form  of  a  flow  diagram  representing  the  computer  model  set  up  in  this  report  is 
shown  in  Fig.  2.  This  carries  out  the  calculation  outlined  in  the  last  paragraph  and  serves  to  introduce 
the  subroutines  MOIN,  MODIEL,  MOMPL,  MOMT,  MOOUT  and  MOPLOT.  We  adopted  a  philoso¬ 
phy  of  breaking  the  main  program  MOFILM  into  a  series  of  subroutines  with  the  knowledge  that 
different  users  would  have  different  computer  facilities  which  might  require  changes  in  MOIN,  MOOUT 
and  MOPLOT.  Also  users  might  wish  to  modify  the  calculation  to  suit  their  own  special  application. 
Such  changes  are  simplified  by  the  modular  approach  used  here. 

The  actual  structure  of  MOFILM  is  somewhat  more  complex  than  indicated  in  Fig.  2.  This  is  due 
to  two  factors.  The  first  is  the  need  for  the  program  to  handle  each  of  the  three  Ken  magneto-optical 
effects.  For  the  transverse  case  MOFILM  uses  the  subroutine  MOMT  when  calculating  the  front  and 
back  surface  scattering  matrices  for  the  magnetic  layer.  Both  the  polar  and  longitudinal  cases,  due  to 
their  formal  similarity,  are  covered  by  MOMPL.  Hence  additional  branching  decision  points  are  needed 
to  distinguish  which  case  is  being  handled  and  to  direct  the  program  to  perform  the  proper  calculation. 

The  second  factor  is  related  to  the  fact,  that  as  shown  in  Fig.  2,  MOFILM  would  make  a  single 
calculation  of  A4>,  R,  \R,  etc.  for  a  fixed  set  of  input  constants.  For  design  purposes,  however,  it  is 
useful  to  see  how  these  (or  other)  parameters  vary  as  a  function  of  changes  in  the  input  parameters. 
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As  a  result,  we  have  set  up  MOFILM  to  permit  easy  variation  of  anyone  of  the  following  parameters:  a) 
the  thickness  of  the  magnetic  layer,  b)  the  thickness  of  the  dielectric  layer  immediately  before  or  after 
the  magnetic  layer,  or  c)  the  angle  of  incidence.  Hence,  suitable  DO  loops  are  superimposed  on  the 
structure  shown  in  Fig.  2. 

Note  that  for  calculations  in  which  the  thickness  of  one  film  is  varied  all  but  one  of  the  scattering 
matrices  remain  the  same.  Hence  when  such  a  case  is  run  the  program  stores  the  product  scattering 
matrices  for  those  parts  of  the  multilayer  structure  before  and  after  the  varied  layer.  This  saves  consid¬ 
erable  running  time  in  determining  the  effect  of  such  a  thickness  variation. 

We  now  discuss  the  main  program  and  the  subroutines  in  greater  detail.  Program  MOFILM  car¬ 
ries  out  the  expanded  version  of  Fig.  2  outlined  above.  Consistent  with  the  ideas  given  earlier,  the 
MOFILM  main  program  is  mainly  a  logical  skeleton  which  directs  the  use  of  suitable  subroutines  to 
make  the  actual  calculations.  Because  of  its  heavily  logical  nature,  particular  care  should  be  exercised 
in  making  alternations  in  the  main  program.  To  aid  in  any  such  alternations,  the  main  program  has 
extensive  COMMENT  statements  to  clearly  define  the  functions  being  carried  out.  At  NRL  MOFILM 
is  run  on  a  PDP-11  computer  in  an  interactive  mode.  Hence  many  of  the  TYPE  statements  in 
MOFILM  are  operator  promptings  and  should  be  treated  as  additional  COMMENT  information.  A 
print-out  of  MOFILM  is  given  in  the  Appendix  as  are  printouts  of  the  subroutines  used. 

MOIN  (NN)  —  This  subroutine  accepts  all  the  needed  input  parameters  and  initializes  the  pro¬ 
gram,  i.e.,  selects  what  M—O  mode  is  being  used  and  what  parameters  (if  any)  are  to  be  varied. 

The  inputs  and  their /unctions  are: 

NO)  —  ambient  index  of  refraction. 

NF  —  number  of  layers  in  structure,  including  the  substrate, 

NF  —  0  used  for  program  escape. 

MO  —  layer  number  of  magnetic  layer. 
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LAMBDA  —  wavelength  in  microns. 

(N,  K)  —  complex  index  of  refraction  n  -  ik  of  the  magnetic  layer. 

(Ci.  Qi)  —  complex  magneto-optical  coefficient  Qx  -  iQi  of  the  magnetic  layer. 

HU)  —  real  index  of  the  Ah  layers,  /  “  1 . NF. 

DU)  —  thickness  of  Ah  layer  in  microns,  /  »  1 . NF.  It  is  not  necessary 

to  insert  thickness  values  for  those  layers  which  have  a 
phase  thickness  of  ir/2  (i.e.,  d  *■  k/An  cos  9) 
since  the  program  automatically  assumes  this.  If  other  fixed 
thicknesses  for  various  layers  are  desired,  they  can  be  entered. 

THETA  —  initial  angle  of  incidence,  in  degrees. 

NTH  —  number  of  THETA  values. 

DTH  —  increment  in  THETA,  degrees. 

L  —  layer  number  of  layer  whose  thickness  is  to  be  varied. 

MTH  —  number  of  L-layer  thicknesses. 

DELMD  —  increment  in  L-layer  thickness,  microns 

NOUT  —  output  control 

NOUT  —  1  printed  output 

NOUT  “  2  graphic  output 

NOUT  -  3  both  outputs 

NOUT  -  0  return  to  beginning  of  input 

NL  TP  —  controls  which  Kerr  effect  used 

NLTP  -  -1  longitudinal 
NLTP  -  0  transverse 


NLTP  -  1  polar 
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NN  —  input  control 

NN  —  0  initial  input 

NN  ■*  1  additional  polarizations 

This  subroutine  also  calculates  the  sine  and  cosine  of  the  angle  of  incidence  in  each  layer  including  the 
complex  values  in  the  magnetic  layer. 

MODIEL(J,DJ)  —  This  subroutine  calculates  the  dielectric-dielectric  scattering  matrix  shown  in  Eq.  (9) 
and  reports  it  back  to  the  main  program. 

J  —  interface  number 

DJ  —  thickness  of  Jth  layer  —  DU),  in  microns. 

MOMPL(M,  NPL,  MD,  D)  —  This  subroutine  calculates  the  magnetic-dielectric  scattering  matrix  for 
the  longitudinal  or  polar  case  and  reports  it  back  to  the  main  program.  The  equations  used  are  shown 
in  Tables  I,  II,  III  and  V. 

NPL  —  control  variable  M  —  magnetic  layer  number 

NPL  »  + 1  polar  MD  —  specifies  front  or  back  surface  calculation 

NPL  *■  -1  logitudinal  MD  —  M- 1  front  surface 

MD  —  M  +  1  back  surface 
D  —  thickness  of  second  layer,  microns. 

MOMT  (M,  MD,  D)  —  This  subroutine  calculates  the  magnetic-dielectric  scattering  matrix  for  the 
transverse  case  and  reports  it  back  to  the  main  program.  The  equations  used  are  shown  in  Table  I,  IV 
and  VI. 

M ,  MD,  D  —  defined  in  MOMPL  earlier 

MOOUT  (NENT,  IN)  —  This  subroutine  calculates  the  multilayer  Fresnel  coefficients  R„,  R„,  R# 
and  Rv  as  well  as  A<b,,  Rh  \Rh  «,  and  the  figure  of  merit  based  on  the  scattering  matrix  for  the  full 
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stack.  MOOUT  also  outputs  the  desired  data  with  the  necessary  labels  in  either  a  printed  or  a  plotted 
format  or  both.  To  perform  plotting  MOOUT  calls  on  MOPLOT.  MOOUT  makes  use  of  the  equations 
in  Table  VII  as  well  as  Eqs.  (10)  to  (12),  Eqs.  (13)  to  (17)  or  Eqs.  (18)-(22)  as  appropriate. 

NENT  —  controls  entry  point  in  MOOUT 

IN  —  keeps  track  of  sign  of  Q. 

MOPLOT  (NLTP,  MTH,  NO,  LM1,  DELMD)  —  This  subroutine  produces  suitable  plots  of  the  vari¬ 
ation  of  output  parameters  such  as  differential  phase  shift,  total  reflectivity,  differential  reflectivity, 
mode  mixing  and  figure  of  merit  as  a  function  of  the  thickness  of  the  magnetic  layer  or  of  nearby  'tun¬ 
ing*  layers.  It  will  also  plot  output  as  a  function  of  the  initial  angle  of  incidence.  Since  this  subroutine 
makes  use  of  local  plotting  subroutines  which  often  vary  from  one  facility  to  another,  it  will  probably 
require  modifications  or  complete  revision  in  order  to  operate  at  another  facility.  The  subroutine  listing 
is  given  in  the  Appendix,  however,  to  illustrate  the  functions  needed. 

NO  —  control  variable  NLTP,  MTH,  DELMD  —  defined  earlier  in  MOIN 


NO  —  2  angle  varied 

LM1  —  control  variable  for  graph  captioning 

NO  —  1  thickness  varied 

LM1  -  0 

layer  just  before  magnetic  layer  varied 

LM1  -  I 

magnetic  layer  varied 

LM1  -  3 

layer  just  after  magnetic  layer  varied 

Certain  mathematical  subroutines  used  in  MOFILM  are  contained  in  our  local  matrix  library  and 
may  not  be  available  at  other  facilities.  We  define  their  functions  here.  It  may  be  desireable  to  substi¬ 
tute  local  variants  for  these. 

CSCLA  (SD,  X,  4,  4,  0)  —  sets  each  element  in  the  complex  4x4  matrix  SD  equal  to  the  complex 
number  X. 

CDCLA  (SD,  X,  4,  0)  —  sets  each  diagonal  element  in  complex  4x4  matrix  SD  equal  to  the  complex 


number  X. 
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CMCPY  (R,  S,  4,  4,  0)  -  carries  out  the  functions  S0  -  Ru  for  each  element  of  the  4  x  4  complex 
matrices  5  and  R  i.e.,  it  copies  R  into  S. 

CMPRD  (R,  S,  T,  4,  4,  0,  0,  4)  —  carries  out  the  matrix  product  T  -  R  x  S  for  the  4  x  4  complex 
matrices  R,  Sand  T. 

This  completes  our  discussion  of  the  details  of  program  MOFILM.  The  printout  of  the  main  pro* 
grams  and  each  of  the  subroutines  MOIN,  MODIEL,  MOMPL,  MOMT,  MOOUT  and  MOPLOT  are 
located  in  the  Appendix.  The  above  discussion  in  conjunction  with  the  COMMENT  statements  and 
some  of  the  TYPE  statements  in  MOFILM  and  its  subroutines  should  make  the  operation  of  the  pro¬ 
gram  fairly  straightforward.  However,  this  chapter  should  be  read  carefully  in  conjunction  with  the 
Appendix  to  avoid  the  common  GIGO  syndrome. 
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IV.  TYPICAL  MULTILAYER  CALCULATIONS 

To  illustrate  the  capabilities  of  program  MOFILM,  in  this  chapter  we  carry  out  a  number  of  calcu¬ 
lations  of  the  magneto-optical  and  optical  characteristics  of  several  magnetic-dielectric  multilayer  thin 
film  structures.  In  particular,  typical  structures  of  the  type  which  might  be  used  as  magnetic  bias  mir¬ 
rors  in  ring  laser  gyros  are  emphasized.  The  results  of  the  calculations  show  the  usefulness  of  the  pro¬ 
gram  as  a  design  tool  for  comparing  various  mirror  designs  or  optimizing  the  parameters  of  a  given 
design. 

Typical  structures  which  might  be  used  for  an  RLG  magnetic  mirror  are: 

A.  An  opaque  magnetic  metal  film  overcoated  with  a  multilayer  dielectric  (MLD)  stack. 

B.  A  transparent  non-metallic  magnetic  film,  such  as  a  garnet,  on  top  of  an  MLD  stack. 

C.  A  partially  transparent  (very  thin)  magnetic  metal  film  on  top  of  an  MLD  stack. 

We  also  discuss  how  the  figure  of  merit  FM,  which  we  define  as 

FM  -  A<&/(1  -  Reflectivity)  (25) 

is  affected  by  the  design  parameters. 

In  order  to  make  things  definite  we  assume  that  the  dielectric  stack  is  made  up  of  layers  of  Si02 
with  a  real  index  of  refraction  n(Si02)  “  1.457  and  ZnS  with  n(ZnS)  —  2.355.  We  ignore  the  small 
changes  in  n(Si02)  and  n(ZnS)  between  A  —  0.63  and  1.15  /xm,  since  the  results  are  rather  insen¬ 
sitive  to  these  changes,  especially  if  the  layers  have  fixed  phase  thicknesses  as  is  true  for  much  of  the 
MLD  stack.  We  also  take  0in„„|  —  30°  and  use  the  values  of  n,  k,  Qh  and  Q2  in  Table  VIII  for  the 
magnetic  layer  under  study.  We  now  discuss  the  above  cases  in  some  detail. 
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Table  VIII.  Parameters  of  Magnetic  Materials  Used  in 
Multilayer  Mirror  Calculations 


Material 

A  (/im) 

n 

a 

Qx 

Qi 

Ref. 

Fe 

1.15 

H 

5.1 

0.053 

E 

“17" 

MnBi 

0.63 

B 

-0.089 

1 

15 

MG3 

0.63 

HI 

0.005 

0.0038 

HSH| 

16 

MG  3  —  (Gd,Pr,Bi)j  (Fe.Ga)j  012  with  —  25,200  deg/cm  and  a  -  4tr k/\ 
estimated  to  be  1000  cm-1. 


A.  MLD  stack  on  top  of  an  opaque  magnetic  film 

The  general  form  we  take  for  this  mirror  is  n(L  -  H)  —  M where  L  is  a  low  index  and  H,  a  high 

index  dielectric  and  AS  is  the  opaque  magnetic  substrate.  All  the  layers  except  that  adjacent  to  AS  are 
taken  to  have  a  phase  thickness  <t>,  —  (2rrn(<f,  cos  0,/A)  —  ir/2  in  order  to  maximize  the  MLD 
reflectivity  for  the  given  We  have  examined  cases  with  an  odd  number  of  dielectric  layers  and 

the  results  are  qualitatively  similar  to  those  shown  here.  In  this  design,  as  the  reflectivity  is  improved 
by  adding  more  layers,  less  and  less  light  reaches  the  magnetic  layer  so  that  the  differential  phase  shift 
A4>  of  the  total  structure  is  decreased. 

Transverse  case  —  Note  that  depends  on  (a)  the  relative  amplitudes  of  the  electric  vector  for 
the  total  light  reflected  from  the  mirror  and  that  part  of  the  electric  vector  of  the  light  reflected  from 
the  magnetic  layer  which  reverses  when  the  magnetization  (and  hence  Q)  reverses  and  b)  the  relative 
phase  between  these  two  vectors.  This  point  has  been  discussed  before  [Report  No.  2,  Chap.  Ill,  Sect. 
A]  [21. 

Macek  [10]  pointed  out  how  one  can  tune  the  relative  phase  by  adjusting  the  thickness  of  the 
dielectric  "tuning"  layer  next  to  the  magnetic  layer.  (We  recall  that  only  a  p- polarized  beam  has  an 
Af—O  effect  for  the  transverse  case  so  all  results  refer  to  the  counterpropagating  p- polarized  beams.)  To 
illustrate  this  effect  we  assume  that  the  mirror  is  6(1  -  H)  -  Fe  with  the  1.15  /*m  Fe  parameters  in 
Table  VIII. 
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The  dependence  of  A4>,  the  p  reflectivity  R,  and  AT?  on  the  tuning  layer  thickness  is  shown 
in  Figs.  3-5.  (The  notation  /I  AT  on  a  scale  means  the  actual  unit  is  10~}  times  that  shown,  for  /£ 5  it  is 
10-5  times,  etc.)  Note  that  there  is  a  "resonant"  behavior  in  all  these  quantities  when  dM-\  =  0.10  /zm. 
Near  this  value  the  first  few  dielectric  layers  adjacent  to  the  magnetic  layer  act  like  an  impedance 
transformer  to  couple  the  light  efficiently  into  the  lossy  magnetic  film.  This  means  that  the  non-Q 
dependent  electric  vector  of  the  multiply-reflected  beam  from  the  magnetic  layer  adds  out  of  phase  with 
that  from  the  adjacent  layers  giving  rise  to  the  dip  in  R.  The  0-dependent  electric  vector  of  the 
multiply-reflected  magnetic-layer  beam  has  a  different  phase  associated  with  it  and  hence  gives  rise  to 
the  enhancement  and  rapid  variation  of  Ad>  and  AT?  for  near  0.1  #im.  Note  that  Ad>  is  maximized 
at  nearly  the  same  <TV_ i  value  for  which  A R  is  zero  as  anticipated  by  the  simple  analysis  [Report  No.  2, 
pp.  9-11]  [21  which  neglects  multiple  reflection  effects.  The  variation  of  the  figure  of  merit  with 
is  shown  in  Fig.  6. 

By  examining  such  results  for  mirrors  with  different  numbers  of  layers  in  the  MLD  stack  one  can 
find  how  (A<t>)m„  varies  with  the  corresponding  total  mirror  loss  (1  -  R)  The  results  are  shown  in 
Fig.  7  where  (A4*)m„  is  clearly  proportional  to  the  mirror  loss.  Thus  (Ad>)m.,/(1  -  R)  is  a  constant  for 
a  given  magnetic  material  (and  9 ,  nL,  nH),  independent  of  the  details  of  the  MLD  design  and  hence 
serves  as  a  useful  figure  of  merit  to  characterize  the  magneto-optical  materials  used  as  opaque  sub¬ 
strates.  This  type  of  analysis  was  used  in  Fig.  9  of  Report  No.  2  [2]. 

Polar  case  —  The  same  type  of  calculations  presented  for  the  transverse  case  were  repeated  for 
the  polar  case  using  the  same  mirror.  In  addition,  the  mode  mixing  eK  [see  Eq.  (19a)]  was  calculated. 
Remember  that  in  the  polar  case  A4>  and  A R  refer  to  the  difference  between  incoming  RCP  and  LCP 
beams.  The  behavior  of  R  is  so  similar  to  Fig.  4  that  it  is  not  plotted.  The  other  results  are  shown  in 
Figs.  8-11.  The  qualitative  features  of  Ad>  and  AT?  as  a  function  of  the  dM- 1  also  resemble  the 
transverse  case.  However,  we  wish  to  draw  attention  to  the  fact  that  one  obtains  a  much  larger 
differential  phase  shift  (and  A R)  in  the  polar  case  even  for  the  same  magnetic  metal,  (cf..  Figs.  3  and  8 
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and  Figs.  5  and  9).  There  is  a  corresponding  increase  in  FM.  The  mode  mixing  eR  in  Fig.  11  is  an 
optical  (not  a  magneto-optical)  effect.  For  the  case  shown  here  the  amplitude  of  the  "wrong"  electric 
vector  is  2.5%  that  of  the  correct  one  at  4*^. 

Finally,  to  illustrate  what  can  be  accomplished  with  an  opaque  polar  Kerr  mirror  based  on  MnBi 
we  show  Figs.  12  and  13.  Note  that  because  MnBi  has  its  magnetization  normal  to  the  plane  of  the  film 
(if  properly  constructed),  it  is  well  suited  for  polar  Kerr  use.  Actually  Fe,  which  was  shown  in  the  pre¬ 
vious  example,  is  not  suitable  for  use  in  the  polar  mode  since  the  shape  anisotropy  causes  its  magnetiza¬ 
tion  to  lie  in  the  plane  of  the  film. 

B.  "Transparent"  magnetic  film  on  top  of  an  MLD  stack 

In  this  case  we  take  the  general  form  of  the  mirror  to  be  M-  n(L  -  //).  Also  we  take  <f>,  -  ir/2 
for  all  the  dielectric  layers  so  as  to  maximize  the  reflectivity.  In  contrast  to  case  A,  in  this  design  all  of 
the  light  impinging  on  the  mirror  reaches  the  magnetic  layer.  As  a  result,  stays  constant  as  one 
adds  more  layers  to  the  MLD  stack.  By  the  designation  "transparent,"  we  mean  that  <<  nM  for 
the  magnetic  film.  In  general,  there  is  some  optical  loss  in  this  film  and,  especially  at  0.63  /im,  this 
may  be  the  limiting  factor  in  the  performance  of  the  mirror. 

The  idea  of  using  combinations  of  "transparent"  magnetic  layers  and  MLD  stacks  as  biasing  ele¬ 
ments  in  ring  laser  gyros  has  been  advanced  earlier  [5,121.  Here  we  merely  show  the  capabilities  of 
MOFILM  for  handing  such  designs  and  illustrate  some  of  the  general  features  of  their  behavior.  For 
the  magnetic  layer  we  take  Bi-doped  garnet  MG3  shown  in  Table  VIII  for  0.63  34m  and.  in  contrast  to 
case  A,  we  vary  the  thickness  of  this  layer  to  properly  phase  match  the  Kerr  reflections  at  its  front  and 
back  surfaces.  Since  it  is  likely  that  magnetic  garnets  will  be  grown  by  liquid  phase  epitaxy  on  a  sub¬ 
strate  of  GGG  (gadolinium  gallium  garnet),  we  assume  that  the  incoming  light  beam  is  traveling  in  a 
medium  having  the  real  index  of  refraction  of  GGG,  viz  n(GGG)  ”  1.95. 
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Transverse  case  —  To  allow  easy  comparison  with  case  A  the  mirror  stack  is  taken  to  be  MG3  - 
6(L  -  H).  Results  for  Ad>,  R,  1R  and  FM  are  shown  in  Figs.  14-17,  respectively.  In  contrast  to  case 
A,  these  quantities  are  slowly  varying  functions  of  the  thickness  varied  ( d m  in  this  case).  The 
differential  phase  shift  for  /^-polarized  light  is  maximized  when  the  effective  optical  thickness  is  nearly 
\/4.  The  slow  decrease  in  R  with  increasing  dM  is  due  to  the  non-trivial  value  of  k  (corresponding  to 
an  optical  absorption  coefficient  a  —  1000  cm-1).  AtX  -  1.15  nm,  magneto-optic  garnets  exist  with  k 
small  enough  to  be  negligible.  The  differential  reflectivity  A R  is  very  small  for  the  case  shown  here 
with  Qi  —  0.  If  one  allows  Q2  to  be  of  the  same  order  of  magnitude  as  Q,,  then  A R  increases  by  a 
large  amount  (see  Fig.  18).  In  fact,  to  a  good  approximation  for  magnetic  materials  with  k  <<  n,  we 
find  Ad>  is  proportional  to  Q\  while  A  R  is  proportional  to  Q2.  Furthermore,  A R  peaks  at  the  same 
value  of  dM  as  does  A<1>  if  Q2  is  significant. 

Although  we  plot  FM  v s  dM  in  Fig.  17,  we  wish  to  emphasize  that  A4>/(1  -  R)  is  not  such  a  good 
figure  of  merit  in  the  case  where  one  has  a  "transparent”  magnetic  overlayer  on  a  MLD  stack.  Since 
Ad>  remains  constant  when  more  layers  are  added  to  the  MLD  while  the  reflectivity  increases, 
Ad>/(1  -  R)  is  a  function  of  the  number  of  layers.  Nevertheless,  one  sees  that  the  phase  shift  per  unit 
loss  is  much  larger  for  the  garnet  mirror  shown  here  than  for  the  Fe-based  mirror  in  case  A.  (If  we 
were  to  compare  both  mirrors  at  0.63  the  difference  would  be  even  more  marked.)  If  k  is  large 
enough  (as  it  is  at  0.63  /xm),  then  A<t>/(1  -  R)  -  FM  does  acquire  a  limiting  value  ’-'hen  more  bilayers 
are  added  to  the  stack  and  its  usefulness  as  a  true  figure  of  merit  is  restored. 

Polar  case  —  The  results  of  calculations  for  a  similar  0.63  /xm  mirror  operating  in  the  polar  mode 
are  given  in  Figs.  19-21.  Again  all  the  calculated  quantities  are  slowly  varying  functions  of  du.  Perhaps 
the  most  surprising  result  is  the  fact  that,  in  contrast  to  case  A,  (A4>)max  is  about  the  same  size  for  both 
the  polar  and  transverse  Kerr  effects  (cf.,  Figs.  14  and  19).  We  point  out  that  in  the  polar  and  longitu¬ 
dinal  Kerr  modes  MOFILM  takes  into  account  the  Faraday  rotation  in  the  magnetic  film. 
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C.  Very  thin  absorbing  magnetic  film  on  an  MLD  stack 

This  case  is  an  extension  of  case  B  to  the  region  where  k  is  comparable  to  n.  This  allows  one  to 
investigate  the  effect  of  Faraday  rotation  and  the  Kerr  reflection  from  the  back  surface  of  a  metallic 
magnetic  film.  The  mirror  structure  is  Fe  -  6  (L  -  H)  and  the  thickness  of  the  magnetic  layer  is  var- 
ried.  The  results  for  this  case  are  given  for  a  transverse  mirror  in  Figs.  22-25.  This  is  merely  illustra¬ 
tive  of  what  can  be  accomplished  with  very  thin  magnetic  films. 

D.  Comparisons  and  comments  on  the  preceeding  cases 

A  typical  40  cm  perimeter  ring  laser  gyro  requires  a  bias  frequency  of  50  kHz  or  more  to  prevent 
lock-in  and  obtain  a  constant  scale  factor.  This  implies  that  the  differential  phase  shift  between  the 
similar  counterpropagating  beams  must  be  >  0.4  mrad  (0.023  deg).  Typically  one  requires  optical 
losses  <  5%  at  1.15  *im  and  <  1-2%  at  0.63  Mm.  The  required  characteristics  of  a  magnetic-dielectric 
mirror  structure  for  RLG  use  are  thus: 

1)  Ad>  >  0.4  mrad 

18  mrad  1.15  pm 

20-40  mrad  0.63  Mm 

3)  No  insuperable  manufacturing  problems. 

Figure  of  merit  values  for  a  number  of  mirror  designs  operating  in  either  the  transverse  or  the  polar 
mode  at  both  0.63  ftin  and  1.15  Mm  are  given  in  Table  IX.  This  table  is  taken  from  the  recent  review 
[3]  of  applications  of  magneto-optics  in  ring  laser  gyroscopes. 

As  yet  no  magnetic  alloy  suitable  for  use  in  a  transverse  Kerr  M—O  mirror  has  been  found  whose 
magneto-optical  performance  significantly  exceeds  that  of  pure  Fe  [3,13],  Table  IX  indicates  that  such 
a  Fe-based  mirror  is  barely  adequate  at  1.15  Min  while  it  would  not  satisfy  requirement  2)  at  0.63  fun. 
On  the  other  hand,  MnBi  operating  in  a  polar  mode  easily  meets  these  requirements  at  0.63  pm. 
Although  we  do  not  have  index  and  magneto-optical  coefficient  data  on  MnBi  at  1.15  Mm.  it  is  likely  to 
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Table  IX.  Magneto-optical  Figure  of  Merit  for  Magnetic 
Metal  and  Garnet  Based  Mirrors 


Mirror  Design 

Mode 

A4>/(opt.  loss) 

X  (/xm) 

MLD  on  Fe 

Trans. 

10 

1.15 

MG4  on  MLD 

Trans. 

>400 

1.15 

MG4  on  MLD 

Polar 

>300 

1.15 

MLD  on  Fe 

Trans. 

6 

0.63 

MLD  on  MnBi 

Polar 

>85 

0.63 

MG3  on  MLD 

Trans. 

230 

0.63 

MG3  on  MLD 

Polar 

210 

0.63 

MG4,  (n  -  2.19,  k  -  1  x  UT\  »F  -  1700  deg/cm)  MG3,  see 
Table  VIII. 


work  well  at  that  wavelength  also.  The  principal  problem  with  MnBi-based  RLG  mirrors  would  seem  to 
be  the  manufacture  of  films  which  actually  meet  the  calculated  performance.  This  is  clear  from  a  litera¬ 
ture  review  of  the  optical  and  magneto-optical  properties  as  well  as  the  physical  structure  of  MnBi  thin 
films  [14,15]. 

One  other  problem  with  the  metallic-based  mirrors  is  the  fact  that  A4>  is  such  a  strongly  varying 
function  of  v  To  achieve  W  =  0  and  Ad>  =  (A4>)m„,  means  that  dM~\  must  be  held  to  tolerances 
of  roughly  50-100  A  if  one  operates  in  the  neighborhood  of  the  "resonance.’  One  can  try  to  get  around 
this  problem  by  operating  away  from  the  "resonance"  with  a  smaller  number  of  layers  in  the  MLD  stack 
to  increase  Ad».  Such  an  approach  allows  Fe  at  1.15  nm  to  remain  marginally  useful  but  does  get 
around  the  thickness  tolerance  problem  for  MnBi-based  mirrors  (see  Figs.  19  and  21). 

The  garnet-based  mirrors  do  not  have  this  problem  because  of  their  slow  variation  of  Ad»  with 
Furthermc  re,  their  figure  of  merit  values  as  indicated  in  Table  IX  are  more  than  adequate.  Their  prin¬ 
cipal  difficulties  are  attaining  an  adequate  Ad>  and  the  fact  that  the  A/l  peak  occurs  at  the  same  du  for 
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which  Ad>  is  maximized.  Although  there  exist  actual  garnets  (for  transverse  Kerr  mode  operation) 
which  have  acceptable  A<b  values,  there  is  little  or  nothing  known  about  their  Q2  values  and  the  conse¬ 
quent  SR' s.  For  the  garnet-based  transverse  Kerr  mirrors  of  Table  IX  the  A<tm„  values  are  0.8  and  2.6 
mrad  at  1.15  and  0.63  fim,  respectively  and  hence  are  adequate. 

The  above  considerations  indicate  that  garnet-based  mirrors  represent  a  good  approach  to  mag¬ 
netic  biasing  for  ring  laser  gyros.  However,  the  good  performance  of  MnBi  in  the  polar  Kerr  mode  at 
0.63  j*m  means  that  it  is  a  viable  alternative  for  constructing  such  mirrors.  For  the  magnetic  garnets  it 
is  desirable  to  develop  materials  with  the  largest  Faraday  rotation  possible  at  the  wavelength  of  interest 
while  keeping  k  and  Q2  as  small  as  possible.  Previous  materials  characterization  work  [16]  has  indi¬ 
cated  that  variants  of  yttrium  iron  garnet  (  YIG,  Y2Fe$0\2)  which  contain  significant  amounts  of  Bi  and 
Pr  plus  Yb  and/or  Gdas  dopants  have  achieved  the  largest  Faraday  rotations  consistent  with  low  k.  It  is 
important  to  characterize  the  Q2  values  of  such  garnets  but  no  data  on  this  exists  to  our  knowledge. 
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In  this  Appendix  we  give  the  computer  listing  of  program  MOFILM  as  well  as  those  of  the  sub¬ 
routines  MOIN,  MODIEL,  MOMPL,  MOMT,  MOOUT  and  MOPLOT.  These  routines  are  supported 
by  certain  matrix  routines  which  are  not  listed  here  but  whose  functions  are  defined  in  Chapter  III. 
Likewise  MOPLOT  is  supported  by  a  number  of  subroutines  which  perform  scaling,  labelling,  tick 
mark,  etc.  functions  involved  in  making  the  actual  plots  of  the  computer  calculations.  It  is  anticipated 
that  local  versions  of  MOPLOT  will  probably  be  written  to  suit  local  variations  in  handling  plots. 

A  study  of  these  listings  together  with  their  functional  descriptions  given  in  Chapter  III  and  the 
COMMENTS  and  TYPE  statements  imbedded  in  the  listings  should  give  the  reader  a  fairly  detailed 
understanding  of  the  overall  program  and  permit  him  to  make  alternations  which  may  better  suit  his 
own  situation.  In  the  final  analysis  all  of  the  useful  optical  and  magneto-optical  information  about  a 
given  magnetic-dielectric  multilayer  structure  is  contained  in  the  complex  components  of  the  scattering 
matrix  (5)  and  hence  any  useful  parameter  describing  the  behavior  of  the  structure  can  be  obtained 
from  ( S ).  This  includes  the  transmission  and  Faraday  rotation  properties  which  are  not  discussed  in 
detail  here  but  which  are  given  in  Table  VII. 

The  various  input  parameters  and  the  calling  parameters  for  the  various  subroutines  are  given  in 
Chapter  III.  Most  of  the  transfer  of  information  among  the  subroutines  and  between  the  subroutines 
and  the  main  program  are  carried  out  via  COMMON  statements. 
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*$•>***  4  v  •>  ■tyrvfiv-t  tv  •>•>  v  v  *  *  *  *  *  v  i  t  ->•/■>*•>  •>•>*<•❖*•>  **♦•>  *  *  *  +  ->  v  ♦  ♦  *  + 


P  !■■:  0  G  R  A  M  H  0  F I L  M 
PURPOSE 

C  A  L  C  U  L  A  T  E  S  M  A  G  N  E  T  0  -  0  P  T I C  A  L  P  R  0  P  E  R  T I E  S  0  F  M  U  L  T I L  A  Y  E  R 
THIN  FILM  SYSTEMS 

S  U  B  R  0  U  "i"  I N  E  S  R  E  0  U I R  E  D 
MGIN(NN) 

MO OUT ' NENT ? IN) 

MO  DIE!...  <  K  7  DK  > 

MOM PL v  M  ? NPL *  MD 7 0 ) 

MOM 7'  (  M  *  MD  7  0  > 

MOP  LOT  1  (  N  •>  M  ?  NO  •>  L  ?  DELX  ) 

THE  CALLING  PARAMETERS 7 EXCEPT  IN  MOPLOT 7  ARE 
NN ? NENT  -  ENTRY  POINTS 

IN  -  BRANCHING  CONTROL 7  RELATED  TO  SIGN  OF  0 
N 7 M 7 MD  -  LAYER  NOS. 

D  K  7  D  -  I...  A  Y  E  R  T  H I C  K  N  E  S  S  E  S 
NPL  -  POLARIZATION 
FOR  MOPLOT 7  THE  PARAMETERS  ARE 
N  -POLAR I Z AT  ION 

M  -  NO,  OF  VALUES  INCREMENTED  VARIABLE 
NO  •-  VARIABLE  INCREMENTED  (D  OR  THETA) 

L  -  LAYER  VARIED 
DELX  -  INCREMENT 

L 1 8 R A R Y  PC U T I N E S  R E Q U I R E D 

COMPLEX  MATRIX  ROUTINES  FROM  DEC  MATH  LIBRARY 
CSCLA 7  CDCLA ? CMPRD 7 CMCPY 
F'  L  0  T  T I N  G  R  0  U  T I N  £  S  F  R  0  H  "f  E  K  T  R  0  N I X 

!. .  0  C  A  L  R  0  U  T I N  E  S 

SCLRND  7  SIZE  7  PAUSE  ARE  USED  IN  MOPLOT 
S  E  E  E  X  P  L  A  N  A  T 1 0  N  T  H  E  R  E 

[  |\t  p  |  J  J  p  £  Q IJ  [  R  £  V'| 

R  E  F  R  A  C  r  I V  E  I N  D I C  E  S  0  F  A  L  L  L  A  Y  E  R  S 

COMPLEX  REFRACTIVE  INDEX  AND  COMPLEX  MAGNETO  OPTIC 
COEFFICIENT  FOR  MAGNETIC  LAYER 


*  4  •»>•>•>  * 


•  *  *  +  *  •>  ♦  4  V  r  r  *  4  ■>  +  •  *  ♦  *  *  »  ■>>»♦■>*  + 

PROGRAM  MOFILM 
IMPLICIT  COMPLEX  (R) 

C 0 M P I... £ X  S P  '■  4  7  4  >  7  D (  4  >  4  )  » S  (  4  7  4  )  *  S 7  M  7  C T M  *  NM  7  Z I  »  D 1 7  <  4  7  4  )  7  g D  <  4  -  4  )  7 1 
1  -  02  T < 4  ?  4 1 

L  A  I...  L  A  M  B  D  A  7  j')  ( 0  )  7  R  K  I...  F1 2  7  R  P  P  P'  2  7  P  P  P  N  2  7  R  P  P  T  2 
R E A L * 8  D D E N T(5 > ?  p 0 L ( 3 ) 

DIMENS  ION  RPP •'  101  )  7 RPS (  101 )  7 RSP  101 )  ?RS3<  101 )  -RRLP2<  101  )  - 


DF'H  I P  1 0 1  >  7  EPS  < :!.  0 1  >  7  DF  EFL  ■  .!.  0 1 )  7  RPP  T  2  <  1 0 .!.  >  7  DP  HI  T  •;  !.  0 1 
2  ,  p  ppp .  1  (>i  )  .  RpF'N  ■;  1  ■;);!  »  ,  $  <  30  >  .  gy  .•  j,j  '■  -  r.  T  <  30  >  7  0  T 1  3  >  ORA  Y 


>  7DELRPU0 
1 12  )  ••  DD  <•  4  > 


'Y>rUN  -  D I.  LAMBDA  7  N  7  0  7  {;•  -  NM  0 !...  T F .  0  •  27  .  S7  ••  C  r'M  ?  STM/SOU T.  S  7  DD  ■  NY  D 
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Vi  :  k  P  L  F/RRLP2  ?  DP  HIP  r  EPS  >  DREFL  »  XO/ ID/ DDE  NT  ,  DT  -  JPEG  T 
2/ I;  I  OUT /'NFO  *  MO  * LO ? NTH ? DTHO t THETAO  r MTH ?  DELMD » NOUT r NLTP 
EQU  !  UALENCE  (  RPP  ( 1 )  r  RPPP 1 )  >  *  ( RI-'S  ( 1 )  r  RPPN  (1))»(RRLP2(1 
1  (  DPH I P  ( 1  >  ?  DPH I T  ( 1  )  )  ?  <  EPS  <  1 )  ■/  DEI..RP  (1 )  > 

D A T A  P 0 L / 'LON G ♦ ' .  ' T RANSU  * ' r ' POLAR ' / 

DATA  DRAY/ ' D ( MA  '  ?  '  G  )  <  , y MU )  ' t ' D < M- '  »  '  1  )  <  r ' MU > ' ? 

I  '  THET  '  r  '  A  (  D  "  ?  EG  >  "  ?  '  .0  <  M+  '  r  '  1  >  (  '  ,  MU  >  '  / 

NN  =  0 

CALL  FOR  INPUT 
CALL  MOIN(NN) 

CON) INUE 
"HE  END? 

I F f  NFO . EQ ,0.)  GO  TO  1 00 

NFIL.~4 

M=MO 

L-L.0 

NF-NFO-1 
DM 1=0 <  L ) 

DM=0 ( M ) 

NNFIL=0 
MNTH=MTH-NTH 
I F  ( MO  » EQ .  0 )  M--NF+2 
LM1  =  1 

if<m.gt.l>  lm  .i=o 

IF(M.LT.L)  LM1=3 

NPllL*NLTP+2 

ODE NT <  5 ) =POL ( NPOL ) 

T  H  E T  A  =  T  H E  T  A  D  * P I ( ) / 180. 

OTH=OTHO*PI ( ) /ISO  * 

CUNT I NUE 

n  ::  wc  ”0 

I F < N N / 2 *2.NE. N N )  N I N C = 1 

NN*NN-NINC 

LOOP  FOR  THETA 

DO  2  I  TH=1 ,« NTH 

ST  a ) ~  S I N ( THETA) 

CT ( 1 >  =COS < THETA > 

CALCULATE  PROPAGATION  ANGLES  AND  THICKNESSES 

DO  3  J=2 r NF+1 

ST  v  J  )  =  ST  <  1 )  *N  1 1 )  /N  <  ..J  > 

C  r  (.  J  >  =SQRT  ( .1 4  -ST  ( J )  M2  > 

IF  •  NN  -  EQ  ■>  2 )  GO  TO  3 
0  <  J ;  -LAMBDA  /  <  A  .  *N  ( J )  *CT  ( J  .■  > 

CONTINUE 

S T M = S T  <  M - 1  >  * N  ( M- 1 )  / N M 
CTM-CSOET  •;  1  *  -S  rM.-X.K2  > 

R  .1  =  <  0  i  • ' 1  *  > 

R  2  ~  J.  *4  A  ■■  1 

IF  f ITH  4  NE » 1 >  GO  TO  5 
CF=1  ^ 

IF  <  N I N  C  ■■■  E  Q .  0  >  C  F  ~  L  AMRO  A  /  ( N <.  L  >  t  C  7  <  L  "•  > 

DM 1 “OMIYCF 
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DELfiD-“CF  4DELMD 
D  <  :l.  :•  -0  * 

.0  (  NFP  :L  )  =0  v 
0M“0 ( M > 

D ( L ) -DM1 

PRINT  OUT  INPUT  INFORMATION 
CALL  MOOUT ( 1 ?0> 

LOOP  FOR  D  ( L. ) 

00  10  JD=1 MTH 

DMaR= ( wD -1 > XDELMD 

I F ( L  . EG ,  M )  Dh“D ( M ) f DMAR 

I F  ■;  L.. ■>  N  E  ♦  1 ♦  A  N  D  ,  L  ■>  N  E  ♦  M )  D  ( L. )  -  0  M 1  4-  D  V  A  R 

0.0  (  l  ■  ~DM 

DO ( 2 ) ~  D ( L ) 

DD < 3 > - 1 80 * * THETA/ PI ( > 

0.0  (  4 )  -0  v  L  > 

IF  uMNTH)  32 ?  32 ?  34 
J-ITH 

£  F  <  I T  H  ,  E  Q ,  1 )  X  0  « 0  0  ( N  X  0  ) 

GO  TO  35 
J-JD 

I F  (  J  0  *  E  Q  ,  1 )  X  0 = 0 .0  <  N  X  0 ) 

I N I ""  I A  L I ZE  M  A  T  R I C  E  S 


call. 

CSCL 

A  ( SO  - 

R1  ?4 

7  4  7  0  ) 

CALL. 

CSC LA ( 3P  ? 

R  .1.  ?  4 

7  4  7  0  > 

CALL 

CSCL 

A  (  D 1 7 

R1 7  4 

«  4  ?  ()  ) 

CALL 

usuL. 

A  <  S  7  R 

j.  ..  4 ... 

4  ?  0 ) 

CALL 

LTJOLA  '■  So  7 

R2  7  4 

7  0  > 

call 

CD  Cl... 

A  <  S  7  F. 

2  7  4  7 

0  > 

CALL 

cdcl 

A  ( SF'  7 

R2  7  4 

7  0  > 

I F  <  NNF I L  .  EQ  ,  0  )  CALL  CSC  LA  <  D 1 T  ••  R 1 . 4  ?  4 , 0 ) 

IF(M  *  EG  *  2)  GO  TO  30 
JFMIN-2 

IF  <  NNF  I L. ,  EQ  -  .1.  )  JFMIN-M-1 
0 1 J  4  J  F  ~  J  F  M I N  M  - 1 
I F  <  N  N  F I L  *  N  E ,  1 )  B  0  T  0  6  6 
C  A  I...  L  C  M  C  P  Y  (  D 1 T  •  S  r  4  *  4  ?  0  > 

DJF-0  ( ,JF  > 

CALCULATE  FRONT  STACK 
C  A  L  L  M  0  D I E  L  ( J  F  ?  D ...  F  ) 

C  A  L..  L.  C  M  P  R  0  (  S  r  0 1 1 S  D  ?  4  ?  4  •  0  ?  0  ?  4  > 

C  A  I...  L  C  M  C  P  V  <  S  0  -  S  ?  4  .-  4  ■  0  > 

IF  (  *  NOT  ,  ( NF IL  ♦  EQ  >  4  .  AND  >  NNF  I L. ,  EQ .  0  ♦  AND .  JF  ♦  EQ  .  <  M-2  ) .  AND  ♦  M » GE  >  4 
1  (30  TO  4 

C  A  L  L  C  M  C  P  Y  ( S  D  ?  D 1 T 1 4  -  4  » 0  > 

NNF  II... -I 
CONTINUE 

IF1  (M— NF)  >  EQ  *  2)  GO  TO  102 
C  H  E  C  K  F  0  R  P  0  i...  A  R I Z  A  T  1 0  N 
IF  < NLTP ;  101 »301 ? 101 

caloiila  -  e  magnetic  layer  for  polar  op  !_ongi  tudinal  polar izat 

C  A  I...  L  M  0  M  P  L  '  M  ?  N  L  I  P  *  M  - 1 1 D  M  ) 
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60 

40 


1. 02 
301 
305 


1 !  0 


CAL L  GAFF'D  (  S D  y  D 1  ?  8  •  4  ?  4  y  0  ?  0  ?  4  ) 

C  A  I..  L  C  M  C  P  Y  (  3  r  S  P  ?  4  y  4  ?  0 ) 

IF  ( ( M-NF )  » EQ *  1)  GO  TO  102 
CALL  MOMPL ( M y  ML TP -M+ly D ( rtf 1 ) ) 

C  A  L  L  CMP  R  D  ( 8  P  y  D 1  ?  S y 4  <4A)rO,4) 

I F  ( . NOT , ’ NF  . GT.M) )  GO  TO  1 02 
IF<NNFIL.NE.2>  GO  TO  105 
CALL  CMCPY <  8  y  SP  y  4  y  4  y  0  ) 

CALL  CMPRD ( SP y  02  T » S y 4 ,4»0» 0 y 4 ) 

GO  TO  102 

C  A  L  L  C  M  C  P  Y  (  S  y  S  P  y  4  y  4  y  0  ) 

CALL  CSCLA ( 8  -  R1 y 4y4r0) 

CALL  CSCLA ( D1 y R1 y 4  y  4  y  0 ) 

C  A  L  L  C  D  C  L  A  ( S  P  R  2  y  4  y  0 ) 

D  0  4 0  J  F  ~  M + 2  y  N  F  + 1 
D  JF--D  <  JF ) 

CALCULATE  BACK  STACK- 
CALL  MOD  I  EL..  (  JF  y  D  JF  ) 

C  A  L.  L  C  M  P  R  0  <  8  y  D 1  y  3  D  y  4  y  4  y  0  y  0  y  4  ) 

CALL  CMCPY ( SO  y  8  y  4  y  4  y  0 ) 

CONTINUE 

C A L L  C M P R D ( SP y 8 0  y  S  y  4  y 4  y  0  y  0 y 4 ) 

IF  (  *  NOT  .  <  NF  I L  ,  EQ  ♦  4  ,  AND  .  NNFIl.  ♦  EQ  ♦  0  >  >  GO  TO  1 02 
C  A  L  L  C  M  C  P  Y  (  S  D  y  0  2  T  y  4  y  4  y  0  ) 

NNFIL-2 

OUTPUT  FOR  POLAR  AMD  LONGITUDINAL  POLARIZATIONS 
CALL  M00UT(2y IN) 

GO  TO  10 
IN~0 

C  A  l.  C  u  L  A  T  E  M  A  G  N  E  T I C  L.  A  Y  E  R  F  0  R  T  R  A  N  8  U  E  R  S  E  P  0  L  A  R I Z  A  T 1 0  N 
CALL  MOM T ( M  y  M- 1 ?  DM ) 

C A L  L  CMP R D  <  3 D  y  D 1  y  8  y  4  y  4  y  0  y  0  *  4  ) 

C  A  L  L  C  M  C  P  Y  (  3  y  8  P  y  4  y  4  y  0  > 

IF  ((M-NF)  <EQ ,  1)  GO  TO  302 
C  A  L  L  M  0  M  T  ( M  y  M  + 1  y  D  ( M  + 1  >  > 

CALL  CMF:'R D  ( SP » D 1  y  S  r  4  y  4  y  0 » 0  y  4 ) 

IF  ( .  NOT  ■>  <  NF  ,GT,  M>>  GO  TO  302 
C  A  L  L  C  M  C  P  Y  (  8  D  y  8  P  y  4  y  4  y  0 ) 

IF < NNFIL . NE » 2 )  GO  TO  110 
CALL  CMCPY ( S y SD » 4  y 4 y 0 ) 

CALL  CMPRD ( SD y D2T y S  y  4  y  4  y  0  y 0  y  4 ) 

GO  TO  46 

CALL  CMCPY ( S  y D2T y  4 , 4  y 0 ) 

CALL  CSCLA ( S - R 1 y  4  y  4  y 0 ) 

CALL  CSCLA ( D 1 y  R1 r 4  y 4  ?  0 ) 

CALL  CDCL..A  (  8  R2  y  4  -  0  ) 

0  0  4  5  J  F  ~  rt + 2  y  N  F  + 1 
DJF  --D  JF ) 

C  A  L  C  U  L  A  T  E  B  A  C  K  S  T  A  C  K  F  0  R  T  R  A  N  S  U  E  R  S  E  P  0  L  A  R I Z  A  T 1 0  N 
CALL  MOD  I  EL  JF  y  D  JF  ) 

C AL L  C rtPRD  ( ,8  y  D 1  y  SD  •  4  y  4  *  0  ?0*  4  ) 

■ :  A  L !...  C  rt  C  P  Y  •:  3  D  y  S  y  4  y  4  y  0  > 
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CON r I NUE 

c  A  1„  l..  C  M  P  E  D  >;  .02 1  ?SD?S?4?4r0r  0  ?  4  ) 

IF  <  *  NOT .  ( NFIL  ♦  EQ  ♦  -4  .  AND  *  NNFIL  ,  EQ  *  0  )  )  00  T  O  46 
CALL  CMCPY ( SD ? D2T  ?  4  ?  4  ?  0 ) 

NNFIL =2 

C A L L  C M C P Y ( S P  *  SD ? 4  ?  4  ? 0 ) 

OUTPUT  FOR  TRANSVERSE  POLARIZATION 
CALL  MOOUT ( 3  *  IN) 

I F  ( I N  ,  EQ .  1  >  GO  TO  305 
CONTINUE 
0  (  L  )  ~  D  H  :L 

IF'MNTH  »LT,  0)  GO  TO  2 

IF ( NOUT  .EQ,  1)  GO  TO  21 

PLOT  -OR  VARIABLE  D 

C A L. L  M 0 P L 0 T  ( N I...  T P ? MTH  ?  1  ?  L M 1  ? D E L M D  5 

CALL  CLEAR 

THETA«THET A+DTH 

IF  < MNTH  » GT »  0)  GO  TO  52 

IF  (NOUT  .EQ.  1)  GO  TO  53. 

PLOT  FOR  VARIABLE  THETA 
C  A  L  L  M  0  P  I...  0  f  <  N  I...  T  P  ?  N  T  H  r  2  ?  0  ?  D  T  H  0  > 

CALL  CLEAR 
NN=1 

MORE  INPUT? 

GO  TO  150 

STOP 

END 


SUBROUTINE  MO IN < NN > 

THIS  SUBROUTINE  ACCEPTS  INPUT 

COMPLEX  Q?NM 

BYTE  YAL 

REAL  LAMBDA ?N 

REALMS  DDENT ( 5 ) 

£■  r  MENS  I  ON  N  <  30 )  -  D  ( 30  >  ••  D  T  ( 3  >  » I  ( 30 ) 

C  0  M  M  0  N  /  D I  /  L  A  M  B  D  A ?  N  ?  Q  -  0  ?  N  M  / 1 D  /  D  D  E  N  T  ?  D  T  ?  J  P  L  0  T 
COMMON  / 1 NOUT/NFO  -  MO  ?  LO ? N  TH  r  DTHO ?  THE  TAD  -  MTH « DELMD  ? 
INOUTrNL TP 

IK(NN.EQ.l)  GO  TO  15 

REENTER  DATA  ?  OR  JUST  CHANGE  OUTPUT? 

CALL  clear 

TYPE  *•'  ALL  INPUT?  EXCEPT  FOR  -ALPHANUMERIC » IS  ACCEPTED •" 
TYPE  BY  LIST  DIRECTED  STATEMENTS' 

TYPE  *?'  TO  ENTER  RUN  I  ,.D,?  TYPE  Y?  IF  NOT?  (CR>  ' 

AC CEP r  152? YAL 
FORMA  <■.  l  A :!.  > 

IF- YAL  . NE .  '  >  QG  ~”0  135 

rYPE  :K?  •  F  'PE  RIJN  IDENTIFICATION  (UP  TO  32  CHARACTERS') 
ACCEPT  42 ?  DDENT 
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FORM  fit!  (  5A8  > 

TYPE  Y?  '  THE  NUMBER  OF  LAYERS  (NF)  IN  THE  GENERAL  CONFIGURATION' 
TYPE  Y  ?  ■'  IS  ARRANGED  AS  FOLLOWS' 

TYPE  Y?'  AMBIENT (1)' 

T  Y  P  E  *  ?  '  D I E  L  £  C  T  R I C  L  A  Y  E  R  S  ' 

T  Y  P  E  *  *  '  M  A  G  N  E  T I C  L  A  Y  E  R  S  ' 

f  Y  P  E  *  ?  '  D I E  L  E  C  T  R I C  L  A  Y  E  R  S  ' 

TYPE  :*?'  AMBIENT  CNF)' 

TYPE  Y?'  ONE  OR  MORE  OF  THESE  MAY  BE  OMITTED ?  EXCEPT  THAT 
TYPE  *»'  THE  TERMINATING  LAYERS  MUST  BE  AMBIENT?  I.E.  ?  OPAQUE' 
TYPE  Y ? '  OR  TRANSPARENT  (NO  LIGHT  SCATTERED  BACK  INTO  SYSTEM)*' 
TYPE  Y?'  NF  «  NO.  OF  FILMS  (NF  =  0  TO  ESCAPE)' 

TYPE  Y " '  M  «  MAGNETIC  FILM  NUMBER  (M  «  0  IF  NO  MAGNETIC  FILM)' 
TYPE  Y- '  LAMBDA  -  WAVELENGTH  (MICRONS) ' 

TYPE  Y ? '  NM  »  COMPLEX  REFRACTIVE  INDEX  OF  MAGNETIC  LAYER' 

TYPE  Y ? '  Q  =  COMPLEX  M-0  COEFFICIENT' 

TYPE  Y ? '  TYPE  NF?  M?  LAMBDA?  <N?K>?  (Ql? 02)' 

A  C  C  E  P  T  Y  ?  N  F  0  ?  M  0  ?  L  A  M  B  D  A  ?  N  M  ?  Q 
I F  ( N  F  0  .  E  0 . 0 .  >  R  E  T  U  R  N 

TYPE  Y? '  NCI)  «  ITH  REFRACTIVE  INDEX  (NOT  ITH  LAYER) ' 

TYPE  Y?'  NN  «  NO.  OF  LAYERS  HAVING  THAT  INDEX' 

T Y P E  Y  ?  '  I (1 ) ♦ . *  I ( N N )  =  L  A Y E R S  H A V I N G  TH A T  1 N D E X ' 

TYPE  Y  ? '  TYPE  NN  ?  N ( I > ? I ( 1 ) ? . * . I ( NN  > ' 

T  Y  P  E  Y  ?  '  (  T  0  E  S  C  APE  L.  0  0  P  ?  T  Y  P  E  /  ( C  R  ) )  ' 

N  ( MO  ) REAL  ( NM ) 

ACCEPT Y  ? I J ? REF ? ( I ( J ) ? J= 1 ? I J ) 

IF(IJ.EQ.O)  GO  TO  11 
DO  5  J=1?IJ 
N ( I ( J ) ) aREF 
I  J  a  0 

GO  TO  12 
CONTINUE 

IF(NN.EQ.O)  GO  TO  10 

TYPE  Y ? '  IF  YOU  WISH  TO  KEEP  PRESENT  VALUES  OF  THICKNESSES 
TYPE  Y?'  TYPE  Y' 

ACCEPT  152  ? YAL 
IFcYAL.EQ. ' Y' )  GO  TO  20 
C  A  L  L  C  L  E  A  R 

TYPE  Y t '  THE  PHASE  THICKNESS  OF  A  LAYER  IS  DEFINED  AS 
T ' PE  * • '  DO ( I ;  a  LAMBDA/ N < I ) YCOS ( THETA ( I >  )  ' 

T  ' !:  ‘  '  (  T  H  E  7  A  ( I )  1 3  T  H  E  F1  R  0  P  Pi  G  A  T 1 0  N  D I R  E  C  T  L  0  N  1  ' 

TYPE  THE  THICKNESS  OF  EACH  LAYER  IS  CALCULATED  AS  DO  k 

T-PE  *? '  IF  YOU  WISH  TO  CHANGE  THE  D'"S?  TYPE  Y' 

TYPE  K  (NO  ENTRY  REQUIRED  FOR  VARIABLE  LAYER  AT  THIS  TIME) 

NN  I.  ~ 0 

0  S  ENTERED  IN  PHASE  THICKNESSES 
ACCEPT  152 > YAL 


IF  •  YAL  .  NE  *  '  Y  '  >  GO  TO  20 

c  * ,  r  >  -TUICKNESS  OF  ITH  LAYER  »  MICRONS  ' 

TYPE  Y  ?  ’  ND  a  NO*  OF  !  AYERS  HA V I NG  THa!  THICKNESS 
T  :  «?  •  I «  1  -■  ,  *  .  I  <  ND  .•  =  LAYERS  Hh1.1  LNG  THA  T  THIGH  NE 

E  k  »  -  TV  PE  ND  ?  D 1  [  '  :•  I  <  1 )  i  t  .  <  [  ■  ND  •  •  CP  •  i  1  'P  EaCH  D  • 
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r  i"F'E  ( 0 '  '  S  NOT  REQUIRED  FOR  AMBIENT  OR  VARIABLE  LA  YE 

L  f-*E  *»'  (TO  ESCAPE  LOOP,  TYPE  /<CR)  >  ' 

AC CEP  I  *» [ J , DTEMP , ( I < J) , J=1 , IJ) 

IF  <IJ.EQ,0)  SO  TO  20 
00  30 v  J=i ,  [J 
O' I (J) >  ~D  TEMP 
[...1=0 
NN 1  =2 

OS  ENTERED  IN  MICRONS 


00 


ro  32 


CALL. 

CLE 

AR 

!  YPE 

K  *  ■" 

NTH  =  NO,  OF  VALUES  OF  THE  I 

A  * 

)  YPE" 

Y  ?  ' 

THETA  --  INITIAL  VALUE  OF  ANGLE 

TYPE 

M  *  ' 

DTH  =  INCREMENTAL  VALUE' 

T  Y  I-1  LI 

Y ,  •" 

L -LAYER  TO  BE  VARIED  IN  THI 

CKNE 

TYPE 

Y  ;  ' 

MTH  =  NO,  OF  VALUES  OF  L  LAYER 

TYPE 

Y  ?  - 

T  Y  p  E  N  T  H  y  T  H  E  T  A ,  D  T  H  ,  1...  y  M  T  H  ' 

AC  CEP 

r  * 

y  NTH .»  THE  TAD  DTHO ,  LO  ,  MTH 

TYPE 

M  ?  •’ 

LAYER  L  IS  USUALLY  VARIED  BY  :!./ 

TYPE 

Y  :•  ' 

DL.0  =  INITIAL  THICKNESS 

OF  i... 

TYPE 

?  f 

01...  1  *  FINAL  THICKNESS  OF 

L. ' 

TYPE 

Y  ,  •' 

NINO  -  0  FOR  DL.0 »  Dl... :!.  IN 

F’  H  A 

i"  YPE 

;X;  „  f 

1 

M .[  C 

TYPE 

Y  ,  '' 

TYPE  NINCy  DL.0  y  DL  L  " 

A  C  C  E:  F 

r  y 

y NINO ,DLO, DLL 

TYPE 

Y y  " 

NOI.JT  =  L  FOR  PRINTED  OUTPUT 

TYPE 

^  ••  •' 

2  FOR  GRAPHIC  OUTPUT 

TYPE 

■M  * 

3  FOR  BOTH' 

TYPE 

M  v  ' 

0  ”0  GO  TO  BEGINNING 

OF 

TYPE 

Y  ,  ' 

NL.TP  =  -1  FOR  LONGITUDINAL 

CASE 

TYPE 

Y  ,  ' 

0  FOR  TRANSVERSE  CA 

SIS' 

TYPE 

Y  ••  ■' 

1  FOR  POLAR  CASE' 

TYPE 

Y ,  ' 

...IP!.,  o  r  =  1  r  0  SUPPRESS  DATE 

ON  G 

TYPE 

? 

0  OTHERWISE" 

TYPE 

If  lt  * 

i-  • 

TYPE  NOUT  ,  NL  IF' ,  JPI..OT  ' 

TYPE 

.'V 

TYPE  0/ < CR  -  fO  RETURN  TO  BEG INN 

Af'CEF 

r  :t 

»NOUT?Nl..  Tpy  JPLOT 

OF  INCIDENCE ' 


SS  ( M  +  :l. ,  M?  OR 
THICKNESSES ' 


PHASE 
US UAL l 


THICKNES 
Y  0 )  ’ 


THICKNESSES 


INPUT 


CAPTION 


OF  INPUT 


D ( LO >  ~0L0 
TE!.MD  =  •:  01...  1  -  01. 
I F  ".  NOU  T  .  EQ 
NN--NN  L  -N  INC 

PS  -"URN 


M  TH- I.  : 
GO  TO 
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S  U  B  R  3  U  T I N  E  M  0  D I E  L. '  K  •  u  K  * 

D I  MENS  I  ON  0 1  (  4  ,  4  >  -  N  <  3 0  >  ,  C  T  (  3  0  - 
IMPLICIT  REAL  ( N  > 

REAL  LAMBDA 

COMPLEX  CX, PROP, PROP  N,D1  t  R 1  , NM » Q ,  CTn  .*  STM 

COMMON  / D I / LAMBDA  ,N»G,D(30),N M /P L  T F / D 1  ,  C T  ,  S  T  (  3 0 )  ,  C T M  ,  S T M 
NC12=N(K-1 )#CT(K) 

NCU=N(K-1  >*CT(K-1 ) 

NC22=N(K)*CT(K) 

NC21*N(K)*CT  <K-1 ) 

PH I ~2 ♦ *P I ( ) * < DK*N ( K ) *CT < K ) /LAMBDA ) 

CX*  <  0  * » 1 ♦ ) TPH I 
PROP~CEXP ( CX > 

PROPN»CEXP  ■;  -CX ) 

RPP*  ( NC21  —NC 1 2 )  /  < NC2.1  +NC 1 2 ) 

RSS* ( NCI 1 -NC22 ) / ( NCI  1 +NC22 > 

TPP-N  (  K- .1 )  *  ( 1 .  +RPP  >  /N  ( K  > 

TSS-1 , +RS3 
R 1  *  ( 0 «  f  0  *  ) 

CALL  CSCLA ( D 1 , R 1 , 4 , 4 , 0 > 

D.t  (  1  ?  1 )  -PROP/TPP 

D 1 ( 1,3) - P R  0  P N  *  R  P P / T P  P 

D1 <  2  7  2 ) ~ P  R 0  P / T  S S 

D 1 (2,4)* P R 0 P N * R S S / T S S 

D.t  (3,1)  *PROP#RPP/TPP 

D 1  (  3  7  3 )  *  P  R  0  P  N  /  T  P  P 

D 1(4, 2  >  *PROP*RSS/TSS 

.01  (4,4)  *PROPN  /  TSS 

RETURN 

END 


SUB R 0 UTI N E  MO M P L  < M , H P L , M D , D ) 

THIS  SUBROUTINE  CALCULATE  SCATTERING  MATRICES  FOR  POLAR  AND 


L  0  N G I T  LJ .0  INAL.  P  0 L  A  R I Z  A  T I C/N  S 

LAYER  MD  IS  BEFORE  OR  AFTER  THE  MAGNETIC  LAYER  IF<M  OR  >M 
RESPECTIVELY 

D I MENS ION  C T ( JO  ) » N < 30 > , ST ( 30 ) 

[ MPL  i.  C I T  COMP!.. EX  <  P .» R  *  ) 

0  !.'  M  E  N  S 1 0  N  R  P  P  ( 2 -  R  S  3  ( 2  > ,  T  P  P  ( 2 )  -  T  S  S  ( 2  >  ,  R  S  P  (  2 )  ,  R  P  S  (  2 ?  T  S  P  (  2 
R E ;  :i L  L  AM9  0  A  ,  N  r  P I 

CO M P L £ X  N C M 1 1  ?  N C M 2 2  ?  N C i 1 1 2  7  N C M 2 1 « N M  7  C T M ,  E M 1  (  4 , 4  )  •  C X  >■  S  f  M  ? 
IDELF'P  ?  DELSS  ?  0 

COMMON  /D 1/  LAMBDA  ,  N ,  0 ,  DU  M  <  30  >  -  Nti/PLTF  -  EMI  ,  CT  »  ST  ,  C  TM  *  STM 
N  C  M :!.  2  -•  N  ( M  D )  *  C  T  M 
NCM 1 :!.  *N  ( MD  >  T C  T  (  MD  ,* 

N  C  M  2  2  *  N  M  >K  C  V  M 
NC  M2 1  «NM-*C  T  ■  M.D  > 

PH  l~2  >  -HP I '  >  KOY.'!'.  MD  >  *C  T  •  MD  •  -  LAMBDA 
l  c .  ^  Q  ,  l  t  ,  rt  )  p  H  L  -  2  -  K  P  I  •  .  V  D  *:  N  M  k  C  T  M  .•  I. .  A  M  B  D  A 
C  •  -  0 , 7  1  .  *PH  I 


y 
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PS  I  -P  I  ( )  :«NM*a*DUh  ( M )  /LAMBDA 
PRP-CEXP ( CX ) 

PRPN--CEXP  ( -  CX  > 

IF(MD-M)  5 . 5  ? 1 0 

Kl~2 

K2-L 

GO  TO  15 
K:L  =  1 

RPP  •;  K1  >  =  t  NC M2.1.  -NCM 1 2 )  /  <  NCM2 1 +NCM 1 2 ) 

E S S(K 1 > ~ ( W C M  Li- N C M 2 2 > / < N C M 1 1 + N C H 2 2 ) 

TPP  (  K .1  >  ~N  ( MB )  M  ( 1 ,  +RPP  ( K  L  >  )  /Or) 

TSS(K1 )=1 .+RSSCK1 ) 

RPP  ( K2 )  -~RPP  ( K 1  > 

RSS( K2 ) “-RSS < K 1 ) 

TPP ( K2 ) -NH* ( 1 » -RPP ( K 1 ) ) /N ( HD > 

TSS  < K2  > -1 » -RSS ( K1 ) 

RSP(K1  )  =  (0,  :.()■>  ) 

I F  ( NPL.  .  EG!  >  - 1  >  GO  TO  1. 

[  F  ( C  T  ( M  D )  .  E  0 .  0  ,  >  G  0  T  0  2 

RSP  <  K 1 )  -  (0.  » .  25  >  *0*  TPP  ( K 1 )  :*TSS  ( K 1  >  *NM/  \  CT  ( HD )  )KN  ( M.D  )  ) 

GO  TO  2 
PSIO=PSI 

P  SI-PSIO* S T h / C  T  M 
I F ( C r ( M D ) * E Q . 0 . >  GO  T 0  2 

R  S  P  <  K  L )  -  <  0  ,»,25)#  Q  %  T  P  P  ( K 1  >  *  T  S  S  ( K 1 )  *  S  T  ( M  D )  /  <  C  T(M  D )  .t  C  T  M  ' 
RPS (  K.'t  )*NPL*RSP(K1 ) 

TSP  (.  K 1  )  -  RSP  ( K 1 ) 

T P S ( K 1 >®(N( M D ) / NM ) # < 1 ♦ - N P L * < 2 , * C T M * C T ( H D > ) / T P P < K 1  )  )  * R P £ 
TAUPP- < TPP < K2 > tTPS < K1 > -RPP ( K2 ) *RPS ( K1 ) ) 

T A U S S - < T S S <  K 2 >*TS P « K 1 ) - R S S 1 K 2  > # R S P ( K 1 ) > 

DELPP=1 .-RPP<K2>*#2 
DELSS- 1 ♦ -RSS  <  K2  >  **2 
IF(MD-M)  20  > 20 ? 25 
PS I 0— PS I 
P  S I  =  C  S I N  (.  P  S 1 0  ) 

R  C  S  -  C  C  0  S  (  P  S 1 0  > 

GO  TO  30 
P  S I  —  <  0  •>  ?  0  *  ) 

RCS:-  ( 1  ,  -  0 .  ) 

R 1  -  ( 0  ♦  » 0  ♦  > 

C  A  L.. !...  C  S  C  L.  A  (  E  H 1  ?  R :!.  ?  4  ?  4  *  0  ) 

EH  1  (  1 1 1 )  -  PCS  -KPRP/TPP  (  2  ' 

Ei '  1  ( 1 » 3 )  - R C S * P R P N R p p  •:  2  >  Tpp  <  2  • 

E  H 1  <  2  •  2  >  -  R  C  S  *  P  R  p  /  T  S  S  ( 2  > 

EH  1  ( 2  ?  4  ) -RCS*PRPN*RSS  ( 2  >  / ""SS  <  2 ) 

EH1  <  3  *  I. )  -  R  C  S  -K P  P  P  *  R  p  p  ( 2  '■  •  T-  p  ;;  ■ 

EH  1.  *  3 » 3  )  - R C S P R P !•! , '' Tpp  ( 2  , 

hTi  L  '  4  ■■  2  )  - R C b M P R P T R S S  '  2  )  ■’  TSS  '  2  1 

EMI  <  4  4  )  =PCS  KPRPH/TSS  (  2  > 

IF' H D — H )  5  3 ■ 3  5  *  •  f  0 

EH  1  <  1 7  2  >  -  PRP  TPS  [ /TPP  ■,  2  > 

TOP-  TPP  <  2  )  fRPS 2  >  -TPS  •  2  )  .HP  pp  (  2  » 


i  <  K 1 ) 
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!"- M 1  ( 1  r  4  )  ■=  (  R C S  ¥  I  UP/  (  TPP  ( 2  >  *TSS  <  2  >  )  —RPP  (  2  ) * P S I  /Jpp  <  p  >  )  F' R F‘ N 

E  M I  ( 2  *  1  >  a  -  P  R  P  #  P  s  I  /  T  S  S  ( 2 ) 
r  i  :i  p  «  r  c  s  #  r  s  p  ( 2  >  +  r  s  s  <  2 )  *  r  p  p  <  2  ;■  *  p  g  1 

E  M  1  (  2 1  3  )  •=  P  R  P  N  *  T  0  P  /  (  T  P  p  <  2  )  *  T  S  S  ( 2  )  > 

EMI  ( 3  1 2 )  =PRP#RPP  ( 2 )  XtPSl  /Tpp  <  2  ) 

TOP  a- ( RCS*TPS  <  2 )  +  TSS  ( 2 )  *F'S  I ) 

EM :!.  (  3  -  4  >  aPRPN*  T  OP/  <  TPP  <  2 )  *T5S  ( 2  >  > 

E  M 1  (  4 » 1 )  -  -  P  R  p  *  R  S  S  ( 2 )  %  P  S  t  /  T  g  g  <  2. ) 

T 0 P = T P P ( 2 > * P 8 1  - T S P ( 2  > :* RCS 

EM  .1.  ( 4  *  3  >  apRPN* TOP  /  ( TPP  <  2 )  * T 8 8  ( 2 )  ) 

SO  TO  45 

EMI  <  1  !•  2 )  «-NPL*PRP*TAUPP/ <  DEL SS# TPP  <  2 )  ) 

E  M 1  <  :l  *  4  )  -  -  N  P  L  :*  P  R  p  N*RSS<2)  t  T  A  U  p  p  /  (  0  £  |__  g  g  #  f  p  f  p  )  ) 

E  M 1(2*1)  =  -  N  P I...  *  P  R  P  *  T  A  ( J  S  3  /  <  D  E  l.  P  *  T  3  S  ( 2  >  ) 

E  M 1  ( 2 1 3  >  =  -  N  p  L  *  P  R  P  N  t  R  P  P  <  2 )  *  T  A  U  S  S  /  <  D  E  L  P  P  *  r  S  S  t  2  >  ) 

E M 1 < 3  f 2 ) ~ NPl.*PRP* R P 8 < 1 } / ( D E L 38* TPP ( 2 )  ) 

EMI  <  3  *  4  >  »NPL  *PRPN*RPS  *  l  )*RSS<2)  /  -  DEL  S3*  TPP  <  2 ) 

EMI  <  4  *  1  >  ajMp|..*pRp*RSp  I  [  ,  /  OELP  P*  I  SS  ••  2  >  ) 

EMI  (  4  ?  3  >  aNPL*PRPN*RPP  <  2  •  *RSP  <  .!.  }  /  <  flp.l  PPtTSS  >  > 

CONTINUE 

RETURN 

D 


S  U  B  R  0  U  r  I N  E  M  0  M  T  <  M » M  D  -  0  ' 

THIS  SUBROUTINE:  CALCULATES  SCATTERING  MATRICES  FOR  TRANS  OF 
POLARIZATION 

0 1 M  E  N  S 1 0  N  C  T  C  3  0 )  » N  ( 3  0  >  .  S  T  <  3  0 ) 

IMPLICIT  COMPLEX  <P*R»T) 

0 1 MENS I ON  RPP ( 2 ) - RSS  < 2 ) * TPP (2) ? TSS < 2 ) 

REAL  LAMBDA  *  N  ?  PI 

COMPLEX  PHI  *  NCM 1 1  *  N  C  M  2  2  ?  N  C  M 1 2  *  N  C  M  2 1  *  N  M  *  C  T  M  *  E  M  :l  {  4*4  >  *CX»BRA 
C 0 M M 0 N  / D I / L A M B D A*N*«* D U M  (  3 0  >  * N M / P L  f F / E M 1  * C T  *ST * C T M  *  S T M 
N  C  M 1 2  a  N  (  M  Ti  >  *  C  T  M 
N  C  M 1 1  a  N  ( M  D  )  *  C  T  <  M  D ) 

NCM22«NM#CTM 
NCM 21  •~-i'lM*C  I  ( MD  ' 

PHI-  2  >  KP I  (  )  *  D  *  N  ( M  D )  >  C  T  <  M  D )  /  L  A  M  B  D  A 

IF  <  MD  .  L  T  .  M  )  PHI  a 2  >  *F  I  (  >  X D t N M >K C T M /  L A M B  A 

C X *.  0  *  9  1  ♦  >*PHI 

PRP-CEXP ( CX ) 

PRPN=CE‘<P  ( -CX ) 

IF  MD  *  GT  ..  M  >  Qa-Q 
B  R  A  a  1 ...  +  <  0  .  1  >  .'  *  Q  *  S  I M  /  C  T  M 
I F  MD— M )  5  ”  '5  ?  1 0 

K  !.  a  2 


GO  TO  .15 


■  •:  J  C  M  2 1  *  B  R  A  -  N  C  M 1 2  )  /  v  N  C  M  2 1  *  &  R  A  4  N  C  n  .1 


RSE 


? STM » 0 
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RSS<  K 1 )  :=  ( N C M 1 1  - N c M 2 2  >  /  ( NCM11+NC M 2 2 ) 

TPP  ( K  1 )  »'N  <  MB  >  *  (  1 .  +RFP  C  K  .1  >  >  / NM 
T5S(K:L  )=i  +  RSS  ( K:L ) 

R P  P  <  K  2  >  =  <  B  R A  *NCM 1 2 - N C M 2 1 ) / ( N C M  2 1 + N C M 1 2 Y ( 
RSS ( K2 ) =-RSS ( K1 > 

TPP  (  K2 )  =NM*  ( 1 .  +  RPP  <  K2 )  ) /N  <  MO  ) 

TSS(K2)»1 .-RSS(Kl) 

R 1  —  ( 0  .  jO» ) 

C  A  L.  L..  C  S  C  L  A  (EMI  r  R 1  r  4 1 4  *  0 ) 

EM  1  ( 1  1  >  ~PRP :*  (  TPP  ( 1 )  #TPP  ( 2  )  -RPP  ( 1  )  *RPP  (  2 
EM 1 ( i ? 3 ) =PRPN*RPP ( 2 ) /TPP 1 2 ) 

E M 1  ( 2  f  2 )  =  P R P / TSS<  2 ) 

EM 1 ( 2 r 4 > =PRPN*RSS < 2 ) /TBS  < 2 > 

EM 1 ( 3 , 1 ) =-PRP*RPP ( 1 ) /TPP ( 2 ) 

EMI ( 3 » 3 ) -PRPN/TPP ( 2  > 

E M  J.  (  4  r  2  >  ~ P RPXcRSS  ( 2 )  / T S S  ( 2 ) 

EM  .1  (  4 1  4  )  aPRPN/TSS  ( 2 ) 

RETURN 

END 


,  ■>  0 ,  )  -BRA  >  ) 


2 ) 


S  U  B  R  0  U  T I N  E  M  (3  0  U  T  ( N  E  N  T  .  I N ) 

THIS  SUBROUTINE  CALCULATES  DESIRED  QUANTITIES  FROM  SCATTERING 
nATRIX  *  ALSO  CONTROLS  PRINT-OUT  * 

IMPLICIT  COMPLEX  (R) 

C  0  M  P  L.  E  X  S  <  4  r  4  >  ?  N  M  -  Z I ,  D 1  T  (  4 ,4>,  0  2  T  (  4  t  4  )  »  S  P  <  4 , 4  >  r  S  D  (4.4)?  0  •  Z 
R E  A L.  L  A M 8 D A » N  (  3 0 )  t  R R L P 2  •>  R P P T 2  ■>  R P P P 2  >  R P P N 2  ?  N M 1  ?  NM2 
R EAL * S  0 D E N T ( 5 ) , P 0 L ( 3 > 

0 1 M  E  N  S 1 0  N  D  <  3  0  >  ?  S  T  (  3  0  >  » C  T  <  3  0  >  »  0  T  <  3 )  *  D  R  A  Y  ( 1 2  > ,  D  D  <  4 ) 

0 1  MENS  I  ON  RPP  ( 1 0 1  )  >  RPS  (  1 0 1 )  ?  RSP  ( .1. 0 1 )  ?  RSS  <  1 0 1  > ,  RRLP2  ( 1 0 1  >  * 

1  D  P  H I  p 1. 0 1 )  f  E  P  S  <  1 0 .1.  >  -  0  R  E  F  I...  ( 1 0 1 )  ?  R  P  P  T  2  <  1 0 1 )  t  D  P  H I T  <  1 0 1  )  t  B  E  L.  R  P  (  1 0 1 ) 

2  •>  RPPP  ( 101  >  f  RPPN  ( 10  .1 ) 

DATA  POL/ ' LONG ♦ "  t  "  TRANSV *' t' POLAR ' / 

DAT  A  BRAY  /  /  0  ( MA  '  ?  •'  G )  (  '  .  '  MU )  '  .  '  B  ( M- '  »  "  1 )  (  '  ?  '  MU )  "  > 

I.  THE  F  *  A  (  0  '  •>  '  EG  )  r  '  D  (  M  f  '  ?  '  1  )  (  '  r  ■'  M U  >  '  / 

C 0 M M 0 N  / D  I / L A M B B A ?M. Q t D t N M / P L T F ' R B U M < 4  * 4 ) ? C T « S T ? C T M » S T M 
COMMON  ''  I NOUT/NFO  MO  ?  LO  »  NTH  r  DTHO  *  THE  TAB  -  MTH  ?  DEL  MB  »  NOUT » NLTP 
L .  APR  PL  f  RRLP2  ?  BPHIP  EPS  *  BRER  L  .  XO/  1 D'  BBEN  )  *  *j  I  ?  JPLOT/SOU  T  ■  S  ••  BD  *  NXB 
EQU I UALENCE  ( RPP  ( 1  >  ?  RPPP  (  I. )  )  ?  <  PPS  v  1 )  -  RPPN  <  l  >  ?  •  \  RRLP2  ..  1  >  ••  RPPT2  (  1  ) 

1  >  -  <  DPH  [  P  ( 1  >  ,  DPH I  T  ( 1  >  )  ■>  ( EPS  ( 1 )  r  DELRP  *  1  ‘  , 

GO  TO  ( 1  *  2 » 3 )  NENT 
PRINT  INPUT  DATA 
CONTINUE 

F 1  (  Z  ’  -  A I M  A  S3  (  Z  )  R  E  A  L  <  Z  ) 

F 2  '  Z  ;  -  AT AN  •  FI  •.  Z)  > 

M~M(. 


L  “  L 1  > 

MM  l  “REAL  ■.  NM  .• 

'  -'-A  [  MhG  '  NM  > 

■••F  “  -'F  )  L 

D  I. .  ” :  '■ :  =  D  ■  !...  •  -  M  T 14  - 1  :•  *  D  E !..  M  D 
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CALL  DATE  IDT)  j 

PRINT  5  •  OT  ?  ODENT  j 

FORMA  r  (  /////  '  DATE  J  '  r  2A4  r  A 1  / 1 X  y  5A8//  )  «j 

PR  INT  30 y  THE TAD  » LAMBDA  ) 

PRINT  35  j 

DO  60  J~ 1 *  NFO  > 

I F  ( J  *  E Q  .  L  .  A N D  .  L.  >  NE  ,  M ;  p R I N T  4 0  ?  J  t  N  ( J  >  t  D  (  J  >  «  D L H  A X 

I F ( J ♦ EQ  »  M  * AND  *  J  *  NE . L )  PRINT  45r J.NM1 *  NM2  r  D <  J  > tQ 

L  F  >.  J ,  EQ  M  »  AND  *  J .  EQ  .  L  >  PRINT  50.  J>NM:L  ?  NM.2  r  0  (  J  >  r  DLMAX  r  Q 

IF  <  J  .  NE  ,  L  ♦  AND  ,  J  .  NE  ,  M  )  PR INT  55  J  •>  N  ( J  >  v  D  (  J  ) 

CONI INUE 

CALCULA TE  PRINTOUT  VARIABLES 

MN7  l-T-MT  H-NTH 

NXDO-O 

IFl.'i.GT.L)  NX  Ei  0=1 
IF-.  M»LT*L)  NX  DO =3 
I F  <  L  *EQ.  1>  NX  DO-- 2 
NXDT=1+3:KNXD0 
N  X  D  T  2 ~ 2 + N  X  D  T 
NXD-l +NXD0 
JAP=0 

RETURN  IF  ONLY  GRAPHICS  REQUESTED 

IF WOUT  ,EQ,2>  RETURN  - 

IF ( ( M-NF >  . EQ ♦  2)  GO  TO  8 
IF'CNL  TP  ..  EQ .  0)  GO  TO  6 

PRINT  J.  000  f  <  DRAY  (  J  )  *  J-NXDT  r  NXDT2  )  * 

PRINT  1100 

RETURN 

PR  INT  2000  >  ( D R A  Y  <  J  >  t  J-NXDT :» NXBT2  > 

RE  TURN 
CENT INUE 
'-■iR  —  JaR  +  I 
,'AR 

AL.  ji.A'fE  Ob  TP'  !T  |;r0R  POLAR  OR  LONGITUDINAL  • '  I  ■',.•••  .A  :D  -  •  • 

-V  \i  :::  C  •  •■.}  ..  I  »  w  <  1  ;  X  ‘  —  £  |.  •>  ;)  »  +! i  4  ?  X  ) 

• J  •  3  -  3  *  r*::.  ■  *  4  •  (  3  ••  4  •  tS  •:  4  ?  *  ;• 

K  p  •'  .  j  f;  ;\l  /  U  [,  ■ 

r ♦  :::  --j  «'  •  «  *  •  3  ?  •:!  '«  *  3  * '3  \  3  «  3  ' 
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I 

I 

1 


I 


I 


! 


W 

i 

i 

i  ■ 


V 


SOS  Pj:  >  5*  (  R1P+R2P  ) 

P  I...  L  0  -  5  *  (  R 1 P  -  R  2  P ) 

R  ft  L  =  0 , 5  *  ( R 1 M  -  R  2  M  > 

R  L  R ™ 0  *  5  #<R1M+R2M) 

GO  TO  310 

308  RRR  -0  ♦  5.*  <  R1P-ZI#R2M  ) 

R R L » 0 , 5  * < R 1M+ZI# R 2 P > 

R  L  R  ~  0 , 5  #  (  R 1 H  -  Z I  *  R  2  P ) 

RLI..=0  . 5*  ( R 1  F'+Z  I  *R2M ) 

3 1 0  PRL2= ( CABS ( RRL ) ) **2 

RLR  2=  <  CABS  (  RLR .»  )  XX2 
RPI...P2  v  J  > a <  0  .  5* CABS  ( R 1 M  >  )  **2 
DPH I P < J )~F2 ( RRL ) -F2 ( RLR > 

EPS  < J ) -CABS < RRR/RLR ) 

DREF'L  (  J  ) » •;  RLR2-RRL2 )  /RRLP2  ( J  > 

I F ( NOUT  . EQ .  2)  RETURN 

PR I NT  1 200 1  DO ( NXB > t RPP < J ) t RPS ( J ) , RSS  < J ) t RSP  < J ) 

PR  I  NT  1 300  f  RRL.P2  <  J  >  ?  DPH  I P  ( J  )  *  EPS  ( J ) » DREF'L  ( J  ) 

RETURN 

3  CONTINUE 

C  CALCULATE  OUTPUT  FOR  TRANSVERSE  POLARIZATION 

3  0 2  R  N  » S  ( 4  y  4  >  *  S  ( t  ?  3 )  -  S  <  1  ?  A  )  X  S  ( 4 , 3  > 

RD=S ( 3 ? 3 > XS < 4 , 4 ) -S ( 3 • 4 > XS < 4 r  3 ) 

1 F  (  ( M  ~  N  F )  .  E  Q  .  1  >  0  ~  -  Q 
I F  -IN  ,  EQ .  1>  GO  TO  303 
JAR- 2  ART  .1 
J5*  JAR 

RPPP < J ) aRN/RD 

IN«1 

RETURN 

305  RPPN  <  J ) aRN/RD 

304  CONTINUE 

I  N:~,.i 

R  P  P  P  2  =  ( C  A  B  S  <  R  P  P  P  <  J  >  )  )  *  *  2 
RPPN 2'"  ( CABS  «  RPPN  (  J )  >  )  **2 
RPPT2  <  J )  *0  5*  ( RPPP2+RPPN2 ) 

DELRP c J  > a < RPPP2-RPPN2 ) /RPPT2 < J ) 

DPH I T < J ) -F2 ( RPPP ( J ) ) -F2 ( RPPN < J ) ) 

IF (NOUT  .EQ.  2)  RETURN 

PR  [NT  2 1 00  r  DD(NXD)  » RPPP  ( J )  r  RPPN  (  J  >  »RPPT2<  J  >  >  OPHIT(  J)  DELPF  <  .1 ' 
"‘ETIJRN 

••  •  >RMA  T  (  '  THETA  ~  '  ?F’6»3?  '  DEG.  '  *  /  '  LAMBDA  «  '  ?F"6  *3*  1  •**  f  •"  <5-.  t.xio 
V5  FORMAT (5X»  LAYER  NO.  '  *3X?  "REFRACTIVE  INDEX  -HXr  THICKNESS 

!.  •.  MICRONS)  •  ) 

4  0  F  0  R  M  A  T '  9  X  ♦  I  2  ?  1 1 X  ?  F  6 . 3  *  1  L  X  r  '  D  M I N  —  '  t  F  7 . 4  *  5  X  y  D  M  A  X  '  F  r  *  4  ) 

>5  FORMA  T '  °  X  ••  [  2 » 3X  *  '  N  a  '  *F6 . 3  *  2X  r  '  K  '  *  F6 . 3  *  1 3X  *  F7 . 4 1 1 6X  *  '  0  ~  '  ? 

I.  2F':> .  4  "■ 

F  i  3  R  M  A  I  (  9  ■ ;  y  1 2  y  3  X  *  "  N  a  '  »F  6  ♦  3  «  2  X  »  '  K  a  ■'  « F  6 . 3  r  3  X  t  "  0  M  I.  N  “  '»F'.  4  * 
t  3X »  UMAX  a  •'  •  F7 . 4  »  3X  ?  Q  ~  '  ?2F9.4) 

'it  "j-  F  0  P  M  I  9  X  >•  1 2  ?  L  1 X  r  F  6.3*  2 1 X  ?  F  '  .  4  ) 

■_  .  j-'  .  ipi-iPT  •  !.  X  .  3a4  y  1 0 X  *  Rpp  '  r  25X  r  '  F>F’S  ?  2-5X  r  RSS  »  25X  ■’  ■  RSF’  1 

cGRMa  r  (24*' .  "R-L  REFL  '  *  I6X?  "DELTA  PHi^RAO  •  -I.7X- 
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!.  EPSILON '  ,2QX?  'DELTA  REFL(L)  '/) 

2 0 0  0  F  0  R  M  A  T  ( /  / 1 X  *  3  A  4 , 9  X  ?  '  R  P  P  P  '  ?  2  4  X  t  '  R  P  P  N  '  ?  1 6  X  r  '  P  R  E  F  L  '  ? 

1 4 X  ?  ■■  DELTA  PHI  ( RAD  )  '  ?2X?  '  DELTA  REEL  '  /  > 

1. 200  FORMAT  ( 1 X  ?  F7 , 5  r  2X  ?  4  <  2E  13.5?  2X  )  > 

1 300  FORMAT ( 20X  ?  4 ( E 1 3 . 5 ? 1 SX > ) 

2  1  00  FORMAT  ( 1 X  ?  F7 . 5  ?  2X  ?  2  <  2E .1 3 . 5  r  2X )  /  3  ( E  1. 2 . 5  ?  2X )  > 

END 


SUBROUT I NE  MOPLOT <  N ?  M  ?  NO  ? L  ?  DELX > 

C  THIS  SUBROUTINE  CONTROLS  THE  GRAPHICS  OUTPUT 

C  ASSUMES  MAXIMUM  OF  1.00  INCREMENTS  OF  VARIABLE  LAYER 

0 1  MENS  I  ON  RRLP2  ( 1 0  I.  >  ?  DP  H I P  <  1 0 1 )  ?  EPSRP  <  1 0 1 )  -  DREFL  ■  1 0  I.  ) 

1 ?  D T ( 3 ) ?  N Y <  9  > , R R ( 4  04)?  C  A P 1 (15)?  X  R  A  Y ( 9 ) 

REALMS  I DENT ( 5 ) ?  YRAY ( 18 ) 

COMMON  /ARRPLT/RRLP2  ?  DPHIP?  EPSRP  ?  DREFL  ?  XMIN/ID/  IDENT .  DT  ?  JPI..OT 
E Q U I V ALEN C E  ( R R  ( 1 )  ?  R R L P 2  (!))?( C A P 1  ( .1.  )  ?  I D E N T  ( 1 )  ) 

DATA  XRAY/  '  D  ( M-  '  ?  M  )  ('  ?  '  MU )  '  ?  '  D  ( MA  '  ?  '  G  )  (  "'  ?  •"  MU )  "  ?  '  THE'!'  ?  '  A  ( D 
1  ?  '  EG )  •"  / 

DATA  YRAY/'PJREFL'  ?  '  ?  'DELTA  PH'  ?  '  I -RAD-  '  ?  "DELTA  RE  ■  'FL'  ? 

1 ' FIGURE  0 ' t'F  MERIT  '  ? ' R  *  L  REEL  '  ?  "  ? ' DELTA  PH  "  ?  'I -RAD- ' ? 

2 '  E  P  S I L  0  N '  ?  "  '  ?  '  D  E  L  T  A  R  E '  ?  '  F  L  '  ?  '  F I G  l.J  R  E  0 '  ?  '  F  M  £  R I T  '  / 

D A T A  N Y / 6? 14? 10? 15? 8? 14? 7? 1 0 ? 1 5 / 

I F ( N  . EQ .  -1)  N=1 
X  M  A X  = 1 . 0 0 * < D E L X t ( M-l ) f  X M I N > 

GO  TO  (5? 10)  NO 
5  D  X  P  2  =  (  X  M  A  X  -  X  H I N  i  / 1 0  . 

C  SCLRND  IS  LOCAL  ROUTINE  FOR  ROUNDING  OFF  NUMBERS 

C  MAY  BE  OMITTED 

.0  X  P  2  -  S  C  L  R  N  D  ( .0  X  P  2 ) 

X  M  A  X  *  X  M I N  + 1 0 .  *  .0  X  P  2 
GO  TO  15 

1 0  D  X  P  2  *  <  X  M  A  X  -  X  M I N )  /  9 . 

DXP 2 “SCLRND < DXP2 > 

XMAX~XM.lN+9  ,  *DXP2 
15  UMAX =5 

IF ( N  . E Q  >  0)  UMAX” 4 

DXF'“0 . 5YDXP2 

DO  20  J=1 ?  UMAX 

KY ~&t  I ABS  <  N  -  f  24  (  J  - 1  >  1 

K  X  ~  6  N  0 + 3  %  L  -  5 

IMIN-t'  J— 1  >*.'!.  01 

IF  (J.EQ.JMAX)  IMIN— 0 

Y  M I N  ~  1 , 

Y  M  A  X - 1 . 

MMAX -M 

I F  •"  N  0  ♦  E  Q .  2  1  MM  A  X  ~  M  -  !. 

IF  ' J ♦ LT  *  UMAX )  GO  TO  45 

DO  40  I®1? MM AX 

RP  f  I )  ~RR  (  I  F 1 0 1  >  /  <  1  ..  -PR  l  ■  ) 

CON  riNUE 
DO  25  I *2? MMAX 
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IF  <  I  *  EQ  » i-iMAX  >  SO  TO  35 
IMINI-XMIN+I 

PREVENTS  SPIKES  CAUSED  BY  IMPRECISION  IN  RATIOS  OF  DIFFERENCES 
IF ( ABS < RR ( I MINI > -RR < IMINI-1)). GT . < 10 . * < ABS ( RR ( IMINI+1 ) -RR < IMINI-1 ) 
1 ) )  )  >  RR  ( I  MIN  I )  =:0 . 5*  ( RR  ( IMINI  +  1 )  +RR  <  IMINI-l )  ) 

Y  M I N = A  M I N 1  <  Y  M I N  ?  R  R  ( I M I N  f  I  )  ) 

YMAX-AMAX 1 < YMAX *  RR < IMIN+I )  ) 

SIZE  IS  LOCAL  ROUTINE  FOR  CHOOSING  SCALES 
NEXT  FOUR  STEPS  MAY  BE  OMITTED 

I F  (  Y  M  A  X  *  N  E  ♦  0 .  >  Y  M  A  X  ™  S 1 G  N  C 1  »  ?  Y  M  A  X  >  *  S I Z  E  <  A  6  S  <  Y  M  A  X  )  ) 

IF  <  YMAX  A..  T*  0.)  YMAX  =0  » 

I F ( Y M I N *  N E*0* >  Y M I N = S I G N < 1 , t Y M I N ) M S I Z E ( A B S ( Y M I N ) ) 

I F  (  Y  M I N  *  0  T  ■>  0  *  )  Y  M I N  ~0  . 

D Y P 2  = <  YMAX- YM I N ) / 1 0 > 

IFC  *NOT\ ( YMIN . EQ . 0 > . OR . YMAX*  EQ  v  0 » >  *  AND  *  YMAX  *  NE * (-YMIN) ) 

1 DYP2--AMAX 1  ( ABS  ( YMAX )  r  ABS  ( YMIN )  >  / 1 0  . 

DYP-0 , 5+DYP2 

C  A  L  L  S  P  L  0  T  S  <  7  *  0  0 1 5  ♦  t  *  7 5 1  *  2 5  .*  X  M  A  X » Y  M  A  X  ?  X  M I N  *  Y  M I N ) 

CALL  BOX ( 0 ) 

CALL  AXIS (2) 

CALL  T I CKS ( DXP ?  DYP > , 1 ) 

C  A  L  L  S  C  A  L  E  (  D  X  P  ? ,  D  Y  P  2  > 

IF  (  JPL..OT  >  48  p  48 » 50 
C  A I... !...  T I T  I..  E  ( C  A  Pi?  4  9  > 

GO  TO  52 

C  A  I...  L  T I T'  I...  E  ( C  A  PI  ,4  0  ) 

0  A  L  L.  X  L  A  B  L  ( X  R  A  Y  ( K  X  >  ?  1 1 ) 

C  A  L  L  Y I...  A  B  L  <  Y  R  A  Y  <  K  Y  )  ,  N  Y  ( 4  )k  N  ~  J )  > 

DO  30  I  1  r  M 
X  a  ( I  - 1 )  %  D  E  L  X + X  M I N 
Y«RR IMIN+I  ) 

I  PL  a  2 

IF  Cl  . EO ♦  1)  £ PL « 3 
C  A  L  L  S  P  L  0  T  C  X  t  Y  ?  I  p  L ) 

CALL  COPY 

PAUSE! >  IS  LOCAL  ROUTINE  FOR  AVOIDING  I NS TRUMENTA!  PROP!  FM 
C  A  L.  L  F’  A IJ  S  E  (10) 

CALL  ALPHA 
CON riNVE 
il ALL  CLEAR 
R ETURN 
END 


43 


KREBS  *ND  MAISCH 


References 

1.  "The  Laser  Gyro,"  F.  Aronowitz,  in  "Laser  Applications,"  Vol.  1,  ed.  M.  Ross  (Academic  Press, 
New  York,  1971)  p.  134. 

2.  "Magneto-Optic  Materials  for  Biasing  Ring  Laser  Gyros,"  "Report  No.  1,"  J.J.  Krebs,  G.A.  Prinz, 
D.W.  Forester  and  W.G.  Maisch,  NRL  Memorandum  Report  3870,  April,  1978;  "Report  No.  2" 
G.A.  Prinz  J.J.  Krebs,  W.G.  Maisch  and  D.W.  Forester,  NRL  Memorandum  Report  4198,  April, 
1980. 

3.  "Applications  of  Magneto-Optics  in  Ring  Laser  Gyroscopes,"  J.J.  Krebs,  W.G.  Maisch,  G.A.  Prinz 
and  D.W.  Forester,  1980  IEEE  Intern.  Magnetics  Conf.,  IEEE  Trans.  Mag.  MAG-16  (1980),  to  be 
published. 

4.  "The  SLIC-7  Laser  Gyro  Inertial  Guidance  System,"  R.F.  Morrison,  E.  Levinson,  and  B.L.  Bryant, 
Jr.,  Proc.  1977  IEEE  Nat.  Aerospace  Electronics  Conf.  (IEEE,  New  York,  1977)  p.  1045. 

5.  "An  Improved  Ring  Laser  Bias  Element,"  R.E.  McClure  and  E.  Vaher,  Proc.  1978  IEEE  Nat. 
Aerospace  Electronics  Conf.  (IEEE,  New  York,  1978)  p.  544. 

6.  "Magneto-Optic  Scattering  from  Thin  Solid  Films,"  R.P.  Hunt,  J.  Appl.  Phys.  38,  1652  (1967). 

7.  "Magneto-Optical  Scattering  from  Multi-layer  Magnetic  and  Dielectric  Films,"  D.O.  Smith,  "Part  I. 
General  Theory,"  Opta  Acta  12,  13  (1965);  "Part  II.  "Applications  of  the  General  Theory."  Opta 
Acta  12,  193  (1965). 

8.  "An  Electrical  Equivalent  Circuit  for  the  Transverse  Magneto-Optic  Effect  in  Thin  Magnetic 
Filins,"  R.E.  McClure,  1980  IEEE  Intern.  Magnetics  Conf.,  IEEE  Trans.  Mag.  MAG-16  (1980),  to 
be  published. 

9.  "Electromagnetic  Theory  of  the  Kerr  and  the  Faraday  Effects  for  Oblique  Incidence,"  C.C.  Robin¬ 
son,  J.  Opt.  Soc.  Amer.  54,  1220  (1964). 

10.  "Ring  Laser  Magnetic  Bias  Mirror  Compensated  for  Non-reciprocal  Loss,"  W.  Macek,  U.S.  Patent 
No.  3,  851,  973  (Dec.  3,  1974). 

11.  "A  Manganese  Bismuth  Magnetic  Mirror  for  the  Raytheon  Laser  Gyroscope,"  T.A.  Dorschner, 
I.W.  Smith  and  H.  Statz,  Proc.  1978  IEEE  Nat.  Aerospace  Electronics  Conf.  (IEEE,  New  York, 
1978)  p.  569. 

12.  "Fabrication  of  Thin  Film  Garnet  Structures  for  Intra-Cavity  Laser  Applications,"  E.C.  Whitcomb 
and  R.D.  Henry,  J.  Appl.  Phys.  49,  1803  (1978). 

13.  "Magneto-optical  Characterization  of  Thin  Films  of  FeX-xBx  ,  Fex_xSix  and  Fe-overcoated  Permal¬ 
loy,"  J.J.  Krebs,  G.A.  Prinz,  D.W.  Forester  and  W.G.  Maisch.  J.  Appl  Phys.  50,  2449  (1979). 

14.  "Magneto-optic  Investigations  of  MnBi  Films,"  G.  Lewicki  and  J.E.  Guisinger,  J.  Appl.  Phys.  44. 
2361  (1973). 

15.  "Determination  of  the  Optical  and  Magneto-Optical  Constants  of  Thin  Film  Manganese  Bismuth 
in  the  Visible  Spectrum,"  R.  Atkinson  and  P.H.  Lissberger,  Int.  J.  Magnetism  6,  227  (1974). 


NRL  MEMORANDUM  REPORT  4364 


16.  "Substituted  Rare  Earth  Garnet  Substrate  Crystals  and  LPE  Films  for  Magneto-optic  Applications," 
M.  Kestigian,  W.R.  Bekebrede  and  A.B.  Smith,  J.  Crystal  Growth  42,  343  (1977).  See  also  other 
papers  on  magnetic  garnets  given  in  Ref.  3. 

17.  The  Q i  and  Q2  values  quoted  here  for  Fe  at  1.15  /tm  are  a  correction  of  those  given  in  Ref.  3. 
Note  that  \Q\  is  not  affected  by  this  correction  and  hence  the  results  given  in  Ref.  3  remain  valid. 


k 


FIGURE  CAPTIONS 

Figure  1  —  A  general  magnetic-dielectric  multilayer  thin  him  structure. 

Figure  2  —  Simplified  flow  diagram  for  program  MOFILM. 

Figure  3  —  Differential  phase  shift  vs  tuning  layer  thickness  d.u- 1  for  the  Fe-based  mirror  design 
shown. 

Figure  4  —  Total  reflectivity  vs  dy-\  for  the  Fe-based  design  shown. 

Figure  5  —  Differential  reflectivity  vs  dy. 1  for  the  Fe-based  design  shown. 

Figure  6  —  Figure  of  merit  A<t>/  (1— R )  vs  4v-i  for  the  Fe-based  design  shown. 

Figure  7  —  Maximum  phase  shift  vs  optical  loss  for  several  Fe-based,  transverse  Kerr  magnetic  mirrors 
at  0.63fim. 

Figure  8  —  Differential  phase  shift  vs  dM„\  for  the  Fe-based  design  shown. 

Figure  9  —  Differential  reflectivity  vs  dw_i  for  the  Fe-based  design  shown. 

Figure  10  —  Figure  of  merit  vs  (  for  the  Fe-based  design  shown. 

Figure  11  —  Mode  mixing  vs  dy- 1  for  the  Fe-based  design  shown. 

Figure  12  —  Differential  phase  shift  vs  dy-\  for  the  MnBi-based  design  shown. 

Figure  13  —  Figure  of  merit  vs  dy-\  for  the  MnBi-based  design  shown. 

Figure  14  —  Differential  phase  shift  vs  the  magnetic  layer  thickness  dy  for  the  garnet-based  design 
shown. 
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Figure  15  —  Total  reflectivity  vs  dM  for  the  garnet-based  design  shown. 

Figure  16  —  Differential  reflectivity  vs  du  for  the  garnet-based  design  shown  (Q2  -  0). 
Figure  17  —  Figure  of  merit  vs  dM  for  the  garnet-based  design  shown. 

Figure  18  —  Differential  reflectivity  vs  du  for  the  garnet-based  design  shown  (Q2  -  Q\). 
Figure  19  —  Differential  phase  shift  vs  dM  for  the  garnet-based  design  shown. 

Figure  20  —  Differential  reflectivity  vs  dy  for  the  garnet-based  design  shown. 

Figure  21  —  Figure  of  merit  vs  dM  for  the  garnet-based  design  shown. 

Figure  22  —  Differential  phase  shift  vs  for  the  Fe-based  design  shown. 

Figure  23  —  Reflectivity  vs  for  the  Fe-based  design  shown. 

Figure  24  —  Differential  reflectivity  for  the  Fe-based  design  shown. 

Figure  25  —  Figure  of  merit  vs  dM  for  the  Fe-based  design  shown. 
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PROGRAM  MOFILM 


Figure  2  —  Simplified  flow  diagram  Tor  program  MOFILM 


6< L-H >  -  FE  J.J5MU  TRANSU 
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—  Differential  reflectivity  vs  dU-\  for  the  Fe-based  design  shown 


6<L~H>  -  FE  1.15MU  TRANSU 


D<H-t>  <MU> 

Figure  6  -  Figure  of  merit  Ad>/(l-/f)  vs  rfw_,  for  the  Fe-based  design  shown 
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Figure  7  —  Maximum  phase  shift  vs  optical  loss  for  several  Fe-based, 
transverse  Kerr  magnetic  mirrors  at  0.63pm 
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based  design  shown 


Figure  9  —  Differential  reflectivity  vs  dw_t  for  the  Fe-based  design  shown 
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6<L-H>  .63MU  02  s  0  TRAHSU. 


Figure  14  —  Differential  phase  shift  vs  the  magnetic  layer  thickness  dM  for  the  garnet-based  design  shown 


6<L-H>  . 63MU  Q2 


Figure  16  —  Differential  reflectivity  vs  du  for  the  garnet-based  design  shown  (Q2  -  0) 


6<L-H>  . 63MU  02-0  TRAHSU. 


Figure  1 7  —  Figure  of  merit  vs  dM  for  the  garnet-based  design  shown 
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Figure  19  —  Differential  phase  shift  vs  dy  for  the  garnet -based  design  shown 


MG3  6<  L-H)  63MU  Q2  «  0  POLAR 


Figure  20  —  Differential  reflectivity  vs  dM  for  the  garnet-based  design  shown 
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Figure  21  —  Figure  or  merit  vs  dM  for  the  garnet-based  design  shown 
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DC  MAG)  (MU) 

Figure  24  —  Differential  reflectivity  for  the  Fe-based  design  shown 


